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Abstract. We review our recent work on linear stability for scalar perturba- 
tions of Kerr spacetimes, that is to say, boundedness and decay properties for 
solutions of the scalar wave equation dgtp = on Kerr exterior backgrounds 
(tM, 3a,J\/)- We begin with the very slowly rotating case \a\ « M, where first 
boundedness and then decay has been shown in rapid developments over the 
last two years, following earlier progress in the Schwarzschild case a = 0. We 
then turn to the general subextremal range \a\ < M, where we give here for 
the first time the essential elements of a proof of definitive decay bounds for 
solutions il). These developments give hope that the problem of the non-linear 
stability of the Kerr family of black holes might soon be addressed. This paper 
accompanies a talk by one of the authors (I.R.) at the 12th Marcel Grossmann 
Meeting, Paris, June 2009. 



1. Introduction 

The stability problem for the Kerr family {M,ga,Ai) of black hole solutions to 
the Einstein vacuum equations 

(1) RM = o. 

the system of nonlinear hyperbolic equations governing general relativity, is one of 
the most important unresolved issues in the theory. 

The ultimate goal would be to understand nothing less than the dynamical sta- 
bility of Kerr as a family of solutions to the Cauchy problem for (1), in the spirit 
of the monumental monograph of Christodoulou-Klainerman [13] (see also [33]) 
proving the non-linear stability of Minkowski space. For the time being, however, 
essentially all work in the black hole case has been confined to the linearised setting. 
The simplest such linear problem is that of scalar perturbations, i.e. the study of 
solutions tp of 

(2) □3V' = 0, 

on a fixed Kerr background {A4,g = ga.Ai)- Equation (2) can be thought of as a poor 
man's substitute for the more complicated problem of gravitational perturbations, 
obtained by linearising (1) around the Kerr family. 

The present paper will review the problem of scalar perturbations (2), in partic- 
ular, the rapid recent progress concerning the Schwarzschild a = and very slowly 
rotating \a\ « M cases. We shall then provide, for the first time here, the essen- 
tial elements of a proof of linear stability for scalar perturbations in the general 
subextremal case \a\ < M, that is to say, boundedness and decay statements (see 
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Theorem 1 and its corollaries) for general solutions of (2) with g = ga,M arising 
from arbitrary finite-energy initial data prescribed on a Cauchy surface. 

In the present introductory section, we shall attempt a self-contained overview, 
which will be fleshed out in Sections 2-7. Sections 8-11 will then give technical 
details of the proof of Theorem 1 in the general \a\ < M case. 

As there is much confusion concerning the state of affairs in this subject, we 
begin with a brief discussion of the classical literature. 

1.1. The classical mode analysis. In the classical analysis, one does not study 
general solutions of (2) but rather individual modes. In the Schwarzschild case 
a = 0, study of this problem was initiated in [38]. A mode is a solution of (2) of the 
form 

(3) <,(r)y„K^,'^)e^"* 

where Ymi denote the standard spherical harmonics, and (i, r, 0, 0) denote standard 
Schwarzschild coordinates. In view of the stationarity and spherical symmetry of 
Schwarzschild, from (2) and the ansatz (3) one derives an ODE for ip'^^iir) in r, 
which can be rewritten in the form 

U + LU u = V u 

where ^ is a nonnegative potential function, depending on m and i. (Remarkably, 
this is also the case for gravitational perturbations [45].) It is immediate that there 
are no solutions of the form (3) with finite energy at t = and Im(aj) > 0. This is 
the statement of 'mode stability' of Schwarzschild. 

The literature on mode solutions (3) for Schwarzschild is now vast and addresses 
questions such as understanding transmission and reflection coefficients (relevant for 
black hole scattering) and locating the complex frequencies w of so-called quasinor- 
mal modes (solutions of type (3) with Im(w) < and finite energy on 'hyperboloidal' 
slices). These have been explored numerically and by methods of asymptotic anal- 
ysis. See [11]. 

The existence of mode- type solutions for (2) on Schwarzschild-the essential start- 
ing point for this classical analysis-follows from the dimensionality of its Lie algebra 
of symmetries generated by the stationary Killing field and the rotations. Re- 
markably, as discovered by Carter, although the Kerr metric for a i= only has 
a 2-dimcnsional algebra of Killing fields, the wave equation (2) on Kerr still ad- 
mits a non-trivial complete separation which allows one again to define mode-type 
solutions of the form 

(4) <,(r)5„,(ac^,0)e""*e-'^* 

so that V'to^C'") again satisfies an ODE that can be rewritten as 

u" + LtP'u - V(aiS)u = 0, 

where the potential V now depends, in addition to the labels (m,^), on the fre- 
quency uj. Here, {t,r,9,(j)) denote Boyer-Lindquist coordinates. The Smi are not 
explict but can be characterized as the eigenvectors of an associated operator de- 
pending on aw, which reduces to the spherical Laplacian if aw = (see Section 8). 
The above separation is related to the complete integrability of geodesic flow, dis- 
covered previously by Carter, and in fact has its origin in the existence of an 
additional "hidden" symmetry of the Kerr metric [44]. This separation property 
was extended to gravitational perturbations by Teukolsky [43] . 
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For the mode analysis of scalar perturbations, the main new feature which 
emerges in passing from Schwarzschild to Kerr is the problem of superradiance, 
first discussed by Zcldovich [49]. The geometric origin of this phenomenon is sim- 
ply that the 'stationary' vector field dt fails to be timelike (see Section 3.1), and 
thus the energy flux associated to this vector field is no longer nonnegative (see 
Section 7.1.1). At the level of mode analysis for w e M, this difficulty reflects itself 
only in the range 

(5) 0<mLo/a< , 

It is precisely in the frequency range (5) that the sign of the energy flux through 
the event horizon of a solution of the form (4) is negative. The transmission and 
reflection coefficients of superradiant modes were probed by Starobinsky [39] , who 
derived bounds in various asymptotic regimes. A priori, superradiance might allow 
for the existence of initially finite-energy solutions of type (4) with Im(a;) > 0, and 
thus, growing exponentially in time. These are excluded (for both the scalar and 
gravitational case) in a remarkable paper of Whiting [47], who thus shows 'mode 
stability' of the Kerr family. Whiting's result can be viewed as the culmination of 
the classical mode analysis of this problem. 

1.2. Linear stability? Despite the successes just outlined, the above classical 
mode analysis was not able to resolve what is perhaps the most basic question 
of linear theory, namely that of linear stability, even in the case of scalar pertur- 
bations (2). Indeed, the statement of 'mode stability' is still completely consistent 
with the statement that general linear perturbations i/j with finite initial energy 
grow unboundedly in time 

(6) limsup'0(t, r, 0, (/)) = oo. 

The reason for the failure of mode analysis to exclude (6) can be summarised 
by the following: 

• It is not a priori clear that solutions of (2) with finite energy can be repre- 
sented as a superposition of modes with w e R. 

• Even if solutions can be so represented, statements at the level of individual 
modes typically do not imply statements for the superposition of infinitely 
many modes. 

Indeed, in trying to 'sum' information at the level of individual modes, one has 
to face certain fundamental physical obstructions which are in fact common to all 
wave equations in more than one spatial dimension and necessarily (sec [30]) lead 
one to consider energy-type quantities. It is only such quantities that will identify 
what it is that bounds ip and excludes (6). In the black hole case, these obstructions 
force us to revisit to what extent we actually understand superradiance, the red- 
shift effect and the role of solutions concentrated near trapped null geodesies, two 
issues which we shall discuss at length later in this paper. As we shall see, it 
is precisely the interaction of these difficulties which gives the present true linear 
stability problem its characteristic flavour. 

1.3. Schwarzschild. A prototype for a true linear stability result is provided by 
the celebrated theorems of Wald [46] and Kay-Wald [27] , which show that solutions 
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ip of (2) on a Schwarzschild background remain bounded by a constant for all time, 
in particular, (6) does not occur: 

Theorem (Kay-Wald [27]). Let ip be a solution of the wave equation on the 
Schwarzschild background {M,g = gi\i) arising from sufficiently regular initial data 
on a Cauchy hyperusrface E. Then 

(7) < C\fD 

in the exterior region up to and including the horizon, where D is a suitable higher 
order energy quantity evaluated on the initial hypersurface S.^ 

The proof of the above Theorem effectively initiated the study of the linear 
stability problem from the modern PDE point of view, according to which: 

• The most fundamental boundcdness statements are at the level of quanti- 
tative control of energy quantities. 

• Pointwisc boundedness statements follow from energy-based statements af- 
ter commuting the wave equation with various vector fields, repeating the 
energy bounds, and applying Sobolcv inequalities. 

In the above, no mode-analysis need be appealed to. Such arguments are essential 
in proving existence theorems for non-linear wave equations, like the Einstein equa- 
tions (1). From the point of view of [46, 27], linear stability away from the horizon 
for (2) essentially thus follows from 'standard arguments' in PDE and the fact that 
the stationary Killing field dt is timelikc in the exterior, allowing the energy method 
to be applied as usual. 

The more novel part of the proof of the above Theorem concerns the horizon: 
Since dt becomes null on the event horizon, its associated conserved energy degener- 
ates there, and this does not allow the standard method to apply to obtain bounds 
on Ip up to and including the horizon (see Sections 3.1 and 5.1). This difficulty was 
circumvented in [46, 27] by a number of clever, but fragile geometric arguments, ex- 
ploiting in a non-trivial way the spherical symmetry but also the discrete symmetry 
of maximally extended Schwarzschild. 

A simpler, and more robust, proof of the Kay-Wald Theorem follows by spe- 
cialising [17], to be discussed below in the context of Kerr, and this new approach 
sheds light on what is the physical origin of the boundedness property at the hori- 
zon: it is precisely the celebrated red-shift effect, which can be quantified by 
exploiting in an appropriate way the generalised energy identities corresponding to 
a suitably-defined local-observer-based energy (see Sections 3.3, 5.3 and 6.2). This 
proof again docs not use any decomposition of -0 into cither spherical harmonics or 
modes. 

An additional linear stability statement supplementing the above Theorem would 
be to prove quantitative decay bounds for solutions, showing not only (7), but that 
Ip can be bounded by a fixed decaying function whose modulus depends only on 
the strength of initial data. (Proving such a statement would be essential from 
the point of view of the question of nonlinear stability, as decay is the only known 
mechanism for nonlinear stability (see the proof of stability of Minkowski space 
[13, 33]).) Even in the Schwarzschild case, such decay results have only become 
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available in the past few years: [8, 14, 16, 34, 31, 42]. Proving these decay re- 
sults has led in particular to better understanding of how various aspects of the 
underlying geometry interact in the analysis of (2), in particular, quantifying the 
'obstruction' to decay caused by the presence of trapped null geodesies associated to 
the photon sphere (see Section 3.2).^ This latter obstruction is effectively captured 
by use of virial-type energy currents whose divergence is non-negative definite and 
degenerates precisely on such geodesies (see Section 5.2). Like the new proof [17] of 
boundedness, the final version of the proof of decay [16] in the Schwarzschild case 
takes place entirely in physical space and does not require frequency decompositions 
or spherical harmonics. (See Sections 5 and 6.) 

1.4. The Kerr family. Up until two years ago, even the analogue of the Kay- 
Wald Theorem was a completely open problem for the Kerr family, in view of 
supcrradiance and the lack of a non-negative definite conserved energy from which 
to start. The first result to appear concerning general solutions -0 of (2), see [17], 
was a proof of the analogue of the Kay-Wald Theorem on a large class of underlying 
metrics a class which included as a special case the very slowly rotating Kerr 
metrics ga,M with jaj « M . Results giving additional decay-type statements, now 
specialised to the exactly Kerr case with jaj « M, rapidly followed in independent 
works (see Chapter 5.3 of [19] and [42]), and later also [3]. We also mention previous 
work concerning non-quantitative statements for azimutlial modes [23, 24]. The 
state of the field up to this point is described in detail in [19]. 

The purpose of the present paper is to review the small |a| « AI proof of [19] 
and complete the picture, by giving, for the first time here, the essential elements 
of the proof of decay results for the general subextremal case |a| < M. (A complete 
presentation, together with a proof of the refinement (13) will be forthcoming in 
[21], [22].) The main resuh is 

Theorem 1. Let ^ he a solution of the wave equation = on the Kerr hack- 
ground g = gM.a 

for \a\ < M with sufficiently regular initial data on a Cauchy 
hypersurface E, and let {Yir}T>o be the foliation defined in Section 4-5. Then, with 
\dip\'^ denoting the square sum of non- degenerate derivatives as in Section 4-6, we 
have 

• Boundedness of energy^ 



(8) 




Integrated local energy decay: for arbitrary r+ < Ri < oo 



(9) 



Jo Je^ 



n{r+<r<Ri} 



J So 



and 



(10) 




Polynomial-time decay of the energy-flux: 



(11) 




This obstruction is not present when one restricts to finitely many spherical harmonics, as is 
often the case in the literature. 

''Note that this energy quantity is thus non-degenerate at the horizon. 
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Here, x is a smooth cutoff Junction which vanishes in a neighborhood of the physical 
space projection of the trapped set, and D denotes the square of an appropriate 
(weighted Sobolev) norm of the data, involving higher derivatives ofip. 

As described above, combining energy boundedness and decay with commutation 
arguments gives pointwise boundedness and decay statements. We give here an 
example of the decay statements that follow: 

Corollary 1. We have the following statement of pointwise decay for ^ 



Before closing this introduction, we give a few comments on some of the ideas 
special to the Kerr case. 

1.5. The interaction of the red-shift, superradiance, and trapped null 
orbits. We have already remarked in the Schwarzschild case the importance of un- 
derstanding the red-shift effect and the role of trapped null geodesies for decay, and 
how these can in fact be captured by energy methods, bypassing completely mode 
analysis, in fact, bypassing any reference to Fourier decomposition and spherical 
harmonics. 

Considering now the Kerr case, there are two main additional difBculties: su- 
perradiance, discussed previously, and the fact that the structure of trapped null 
geodesies is more complicated, at least when viewed purely in physical space. These 
two problems are in some sense coupled, and this coupling can be viewed as an ad- 
ditional, third difficulty. (See Section 7.1.) It turns out that to understand these 
difficulties, it is indeed useful to reconsider Carter's formal separation of ijj into 
modes, but viewed from a slightly different perspective. 

Let us note that the first issue mentioned in Section 1.2 (associated to whether 
general solutions ip can be written as superpositions of modes) can be sidestepped 
with the help of a cutoff function, by restricting consideration to a finite time slab. 
In such a slab, one can indeed represent a given finite-energy solution ip of (2) as 
a superposition of modes (4) localised at fixed frequency-triple w 6 R, m, I. If one 
is able to recover a suitable quantitative estimate on the solution in this timeslab, 
then by a bootstrap argument the estimate can easily be extended to be valid for 
all time. 

Recall from our brief discussion of the Schwarzschild case that the key to showing 
quantitative bounds is constructing suitable energy currents which capture the usual 
energy conservation, the red-shift, and the obstruction of trapped null orbits. Key 
to the original boundedness result [17] on Kerr were two observations: 

• In the case \a\ « M , superradiant frequencies are not trapped. 

• When superradiance is controlled by a small parameter, then it can be 
absorbed with the help of the redshift. 

The two observations taken together allowed (see Section 7.2.1) an understanding 
of boundedness without understanding trapping. In fact, the boundedness results 
were obtained not just for exactly Kerr metrics but for a wide class of spacetimes 
whose metric is C^-close to Schwarzschild. The decay result [19] (which is restricted 
to Kerr spacetimes) in the |a| « M case replaces the first observation above with 




[r^l < C^/DT-\ 



\rip\ < C\fDT-^. 



iV']<a,\/:DT-i^''. 
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a complete treatment of trapping. Here the full potential of the microlocalisation 
provided by (4) is used: In the range of frequencies relevant for trapping, virial 
currents arc chosen separately for each frequency-triple (w, to, £), so as to degenerate 
at a unique value of r depending on this triple. This r- value can be related to a null 
orbit with conserved quantities which arc in turn associated to the frequency triple. 
Behind this construction is the close connection between the separation of the wave 
equation and the separation of the Hamilton- Jacobi equations. (See Section 7.2.2.) 

The main new idea necessary for the general |a| < M case is to revisit also the 
first observation above. In the small \a\ « M case, the fact that superradiant 
frequencies are not trapped follows essentially from the fact that all future-trapped 
null geodesies eventually leave the ergoregion. For general \a\ < M the latter is no 
longer the case. Nonetheless, remarkably, 

• Superradiant frequencies arc not trapped for the entire \a\ < M range! 

In the special (marginal) case of frequency to = 0, the above is related to the fact 
that there are no trapped null geodesies orthogonal to dt- This latter fact plays a 
role in recent work on black hole uniqueness [1]. 

The above observation allows us to adapt the second observation from before to 
the case where superradiance is not a small parameter and still couple the conserved 
energy estimate, the redshift and understanding of trapping (the latter again via the 
microlocalisation achieved by the classical separation) so as to obtain the desired 
estimates (see Section 7.3). 

The resolution of the linear stability problem in the \a\ < M case thus brings us 
full circle. Though based on energy methods, the proof reconnects with the classical 
mode analysis, now considered from a more sophisticated point of view as a tool for 
the 'microlocalisation' of energy currents. Moreover, the method highlights certain 
properties of individual modes (4), indeed of the potential functions V{uj), which 
do not appear to have been noticed previously, but turn out to play a fundamental 
part in the stability mechanism. Thus, we hope that this argument offers something 
new even for the reader impatient with some of the more technical difficulties of 
the proof. 



1.6. Outlook. With Theorem 1 and its corollaries, the problem of linear stabil- 
ity, at least for scalar perturbations, has been essentially completely understood, 
leaving only the extremal \a\ = M case. In view of recent work of Aretakis [4] on 
extremal Reissner-Nordstrom, for extremal Kerr one in fact expects instabilities on 
the horizon. Wc should mention that electromagnetic perturbations on Schwarz- 
schild have been treated in [5] , while the problem of gravitational perturbations is 
the subject of recent work of Holzegcl [26]. 



2. The Kerr space-time 

The detailed description of the maximally extended Kerr space-time {A4,gM^a) 
can be found for instance in [9], [19]. Here we simply recall the local coordinate 
expression of the metric in the regions which will be of interest to us. 
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2.1. Boyer-Lindquist coordinates. The Kerr metric in tlie domain of outer com- 
munication TZ c A4 (corresponding to one of the ends) is given in so-called Boyer- 
Lindquist coordinates {t,r,6,(l)) eRx (r+,oo) x [0,7r] x [0,27r) by 

(14) OM a = -^(dt - asin^ ed(l)f + ^dr^ + p^dO^ + (adt - (r^ + a'^)d<j>f . 

Here 

A = r^-2Mr + a^, = + cos^ 6*, 

and r+ is the largest root of the equation A = 0, 



r± = M ± \/M2-a2; 

The coordinate system restricted to 6* =;t 0,7r, (j> + defines an almost global chart 
on TZ. The coordinate vectorfields T = dt and <i> = 9^ extend to globally defined 
Killing fields and characterize Kerr as a stationary, axi- symmetric space-time. 
Setting a = 0, we recognize the Schwarzschild metric 

(15) gM = -(l-^)d<' + (l-^) \r^ + r\de^ + sin^6d(l)''), 

in standard form. 

It is useful to define a rescaled version of the r coordinate on r > r+ by 

after suitably chosing a value for r* = 0. We will sometimes consider the {t, r* ,9, 0) 
coordinate system. Note that these now range in M x R x [0, tt] x [0, 2tt). 

2.2. Kerr-star coordinates. For \a\ < M, recall that the boundary dTZ in the 
maximally extended Kerr spacetime is a bifurcate null hypersurface Ji = T-L^^Tir , 
where 'W' are known as the future (rsp. past) event horizons. 

The geometry of T-L is not easily understood in Boyer-Lindquist coordinates, 
which degenerate as one approaches this hypersurface. Defining, however, so-called 
Kerr star coordinates (r*,r, 0,0*) on TZ. by 

t*=t + t{r), (/)* = (/) + 4>{r) 

where 

^(r) = (r2 + a^)/A^ ^(r) = a/A, 
ar ar 

then these coordinates extend regularly to be defined on a larger subset of Ai, 

covering in particular "H"*" ^ modulo the standard degeneration of the spherical 

coordinates on S^. In these Kerr star coordinates, the set \ T-T now simply 

corresponds to the set {r = r+}. We note that the Killing fields again take the form 

3. Geometric phenomena of Schwarzschild and Kerr 

In this section we briefly review the geometric features of Schwarzschild and 
the Kerr space-time which will be particularly important for understanding the 
behaviour of solutions the wave equation. 



THE BLACK HOLE STABILITY PROBLEM FOR LINEAR, SCALAR PERTURBATIONS 9 



3.1. The stationary Killing field and the ergoregion. We have aheady re- 
marked that for the Kerr family, the Boyer-Liiidquist coordinate vector field T = dt 
is Killing. We consider briefly its geometry. 

We begin with the Schwarzschild case. From (15), we immediately read off 
that T = dt is timelike everywhere in the domain of outer communications TZ. On 
"H^ u H", T becomes null, vanishing in fact on H'^ n "H". (In the black hole region, 
T becomes spacelike; this is of no concern to us here.) 




As we shall see, the fact that T becomes null is already a difficulty when we 
discuss the notion of energy in the sections to come. In the case of the Kerr 
metric gM,a for a + 0, the situation is far worse! Now, the vectorfield T fails 
to be everywhere timelike in the domain of outer communication TZ. The subset 
£ c int(7?,) where T is spacelike is known as the ergoregion. The boundary d£ of £ 
(in the topology of TZ) is called the ergosphere. 




Recall also the axisymmetric Killing field $. We note that the combination 

L=T+ ,° , ^ 

is null and normal to Ti^ \ Ti^ . Thus Ti^ \ Ti^ is a Killing horizon. The Killing field 
L is sometimes known as the Hawking vector field. We note the identity 

(17) VlL=kL, 

where 

The quantity k is known as the surface gravity. We note that in the extremal case 
\a\ = M (not considered here!), k vanishes. 
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3.2. Separability of geodesic flow and trapped null geodesies. The high 
frequency behavior of solutions to wave equations is intimately related to the prop- 
erties of null geodesic flow. We thus consider these properties in the black hole 
case. 

Again, we begin with the Schwarzschild case a = 0. In view of the dimensionality 
of the span of the Killing fields T, ili, £72, i^s, it is immediate that the Hamilton- 
Jacobi equations separate, and thus, geodesic flow can be understood completely 
(cf. the discussion of separation of the wave equation in Section 1.1). Geodesic 
fiow in the Schwarzschild metric is described in detail in many textbooks. One sees 
from the resulting equations that there are null geodesies 7 which for all afRne time 
remain on the hypersurface r = 3A/, the so-called photon sphere. 




With reference to suitable asymptotic notions defining 'future null infinity' I^, we 
can make the following more general statement: If 7(5) is a future-directed inex- 
tendible null geodesic in Schwarzschild (7?,,(7a/.o) with 7(5) e 72. \ for all s > sq, 
such that moreover for all p e I^. 3^ such that 7(s) ^ J^{p), then lims^oo = 
3M. 

Turning now to the general Kerr case, remarkably, as discovered by Carter [9], 
geodesic flow admits, besides the conserved quantities associated to the Killing fields 
T and <&, a third non-trivial conserved quantity (the 'Carter constant'). Thus, 
geodesic flow remains separable and can be completely understood. Again, this 
separation is related to the separation of the wave equation, discussed in Section 1.1. 
The dynamics is described in some detail in [11]. 

In contrast to the Schwarzschild case, there are now null geodesies with constant 
r for an open range of Boycr-Lindquist-?' values. However, restricting to geodesies 
sharing a fixed triple of the nontrivial conserved quantities, there is at most one 
Boyer-Lindquist r-value, depending on the triple, to which all null geodesies neither 
crossing H^, nor approaching , must necessarily asymptote to towards the future. 

3.3. Red-shift effect. An important stabilizing mechanism for the behaviour of 
waves near black hole event horizons is what we shall here call the "horizon- 
localized" red-shift effect. 

Recall that this is the red-shift relating two observers A and B, where B = Vt(^) 
for r > 0, both crossing the event horizon. In the geometric optics approximation, at 
horizon crossing time, the frequency of waves received by B (measured with respect 
to proper time) from A are damped (in comparison to the frequency measured by 
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A) by a factor exponentially decaying in r. 




The above effect depends in fact only on the positivity of the surface gravity, as 
was established for non-extremal Kerr in Section 3.1. In the extremal case, we note 
that the above horizon-crossing time red-shift degenerates. 

4. The energy method 

In this section we introduce the machinery of compatible currents. This is a 
robust and invariant approach to generating L^-based identities, which can be ex- 
ploited to establish links between the geometry of space-time and the behaviour of 
linear (and often even nonlinear) waves. 

4.1. Compatible currents in physical space. Given a Lorentzian metric g, let 
^ be sufficiently regular and satisfy 

□s* = F- 

We define 

Given a vector field and a function w, we will define the currents 
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Applying the divergence identity between two homologous spacelike hypersur- 
faces 5*1, ^2, bounding a region S, with 5*2 in the future of 5i, we obtain 

(19) _^^j;:[*K^ + J^K^[*] + £:^[*] = _^^Jr[*K.> 
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where ng. denotes the future directed timehke unit normal. A similar identity holds 
replacing V by V,w. 

4.2. Noether's theorem. The most celebrated use of (19) is in the case where V 
is Killing, and = -0 where i/' satisfies Hg'ip = 0. One obtains then the conservation 
law 

(20) ^^j;:m< = l^K^i^Kr 

This is a version of Noether's theorem. 

4.3. The dominant energy condition. When V is timelike, then one has more- 
over that Jj|^[^']'^5. is nonnegativc definite, in fact 

(21) j;:[*K>c(^)Eia„*p 

Q=0 

where da denote derivatives with respect to an arbitrary system of local regular 
coordinates. Thus, given a timelike Killing field, (20) and (21) yield an a priori 
estimate for the norm on 5*2 of all local derivatives of solutions = ip oi the 
wave equation, in terms of a similar quantity evaluated on S'l. 

One should not think, however, that this is the only use of identity (19). If V is 
timelike, but not Killing, then one can estimate the bulk term 

K^[*]<c(5)Pj;:[vI/]n'j 

in terms of the boundary. For solutions of the wave equation, this allows one to 
prove that the left hand side of (19), integrated over S2{t) grows at most expo- 
nentially in an appropriate notion of time t. This is fundamental in proving local 
existence type results. 

4.4. Virial identities. Fundamental for the considerations of the present paper 
will be yet another way of viewing identity (19). While in Noether's theorem and 
in the proof of local existence, one often considers the bulk term K^[5'] as an error 
term, (19) is also very useful when one can show positivity properties for the bulk 
term K^[4']. Then, provided that the boundary terms are also controlled, (19) 
can be used to show the boundedness of a spacetime integral quantity. This can 
be thought of as a weak statement of decay. This is reminiscent of the use of the 
virial theorem. 

4.5. The hypersurfaces Et-. Returning now to the Kerr geometry, we shall typi- 
cally apply (19) in regions bound between the event horizon, null infinity, and the 
leaves of a foliation associated to hypersurfaces J^r to be defined below. 

Let r+ < r„ < i?„ < oo. The hypersurface Sq is constructed as follows. Define the 
function / = f{r*) according to 

firl= f 
Jo 

Then 

So = {t + /(r*) =0}, for r+<r<r, 

So = {i = -/«)}, for r„<r<i?„, 

So = = fin - (fir:) + fiK))}, for Rn < r. 



APA 
(r^ + a^)^ 
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Now let denote the 1-parameter family of diffeomorphisms generated by the 
Killing field T and define the hypersurfaces Sr '■= Vt(So). (Equivalently, we replace 
everywhere t by t - r in the definition of Sq above.) 




We may now denote by t the function on J^(So) which is constant on the level 
sets YiT- The hypersurfaces Sr are space- like: 

ff(VT,VT)< 

Near null infinity 

T = t-r* + 0{M'^r''^), r^oo 

while near the horizon 

T = t* + 0{M^), r^r^. 

By the coarca formula, wc easily see that, denoting by B the region bounded by 
Sti and T,T2 and the ^(^^ ^.^y ^{ri T2)' have that for an arbitrary function F, 



When the measure of integration is not denoted, the volume form is always meant. 
We will use this relation without comment in what follows. 



4.6. Non-degenerate derivatives. Finally, since energy estimates will control 
various combinations of derivatives, it is useful to have suggestive notations for 
these. 

First, we shall denote ^ the induced gradient on the r =constant spheres. 

We shall use the notation I^VP to denote the type of quantity controlled by the 
flux of through T,r for r > 0, where is a (^(-invariant vector field which is 
strictly timelike and N = T for large r (for instance, we may take the vector field 
of Proposition 5.3.f ). Explicitly, wc may define 

\d^f = {dr^Pf + {dt*^f + \nf 

for r < r„, where the coordinate derivatives arc interpreted in Kcrr-star coordinates 
and 

\di:f = a{r){dt*ijf + {dr^f + m^ 
for r > r„, where a(r) is a positive function which vanishes as r ^ 00. 
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Finally, let us use the notation \d%j:\^ to denote the derivatives appearing in the 
flux of an N as above on the null hypersurfaces "H^ and I^, respectively. These are 
all tangential derivatives. We may write explicitly 

= \LM^ + \nt 

on T-L^ and 
on I+. 

5. SCHWARZSCHILD: THE THREE SOURCES AND COMPONENTS OF DECAY 

In this section, we describe three ingredients of the recipe for the robust deriva- 
tion of the quantitative decay rates for solutions of the wave equation on Schwarz- 
schild black hole space-times. These are directly related to the geometric phenom- 
ena studied in Section 3. We will see in the following section how these ingredients 
are used to show 'integrated local energy decay' and, together with an extra ingre- 
dient of a much more general type, actual decay. 

The Kerr problem to be addressed later will require revisiting the strategy and, 
among other things, microlocalising the components of this section in a manner 
strongly attached to the Kerr geometry. See Section 7. 

5.1. The conserved energy estimate. The Schwarzschild metric is stationary 
so we may apply Noether's theorem (20) for V = T, where T is the 'stationary' 
Killing field. Specifically, we apply (20) to the space-time region bounded by two 
space- like hypersurfaces St-^ , E^-^ , part of the event horizon 'H^^^ ^ and null infinity 
I^L ^ -,, i.e. we take 52 = "Ht ^ x u E^, ult ^ 5i = Et-, . We obtain: 

(22) ^Imn'^^^J^^ J^[*K, + ^^ Jn*K. = X J^[*K.- 

We shall see in Section 7.1.1 that for the general Kerr case, the above identity 
is a priori of limited use in view of the failure of T to be timelike in the exterior. 
Already, however, in the case of the Schwarzschild solution with a = 0, the vectorfield 
T, though timelike in the domain of outer communication 7?., becomes null on the 
horizon Ti.^ = {r = 2M}. Thus, while all integrands in (22) are nonnegative and 
thus each integral taken individually is controlled by the right hand side, the first 
two integrands degenerate with respect to what one would have were T timelike. 
For the first integrand, we thus only have the degenerate lower bound 



(23) 3lmn'^^{r)>c{r)m\ 

where c(r) = for r = 2M , and similarly, for the second integrand, we only have 

(24) Jn*]"«- ^ 

i.e. the right hand side is missing"' the term j^K'p. 
From (22), we thus obtain the estimates 

(25) f c{r)\dM'< f JH*]^"' 

(26) Li^^i'^ L jn*]"^ 



^cf. the 'expected' 



THE BLACK HOLE STABILITY PROBLEM FOR LINEAR, SCALAR PERTURBATIONS 15 



(27) f \dyj\'< f J^[*]r.^ 

The standard approach to proving pointwise boundedness of solutions ■0 of (2) is 
by proving L^-estimates for all derivatives of ip up to sufficiently high order. Thus, 
even in the Schwarzschild case, the above degeneration of (25) or (26) presents 
significant comphcations for proving stabihty at the horizon. See Section 6.2. 



5.2. Virial estimates and trapping. For the discussion to follow here, it is useful 
to compare with Minkowski space. Defining a current J'''"'"' [4"] where 

(28) X = f{r)dr 

for specific choices of /(r) and u'(r), one has that for suitable regions B 

(29) f K^"-"'[*]>Cj? f (|a*|2 + |*|2), 

JB JBn{r<R} 

and thus, for solutions ■0 of the wave equation, identity (19) yields 

(30) f , m\' + m<C f J^^-Mn^ +C f Jf-Wn^ . 

JBn{r<R} J Si ^ J S2 

The above relation together with energy conservation (see Section 6.1) will al- 
ready yield a weak statement of decay. This for instance already excludes stationary 
solutions. 

It turns out that for more general metrics, trapped null geodesies are a "high 
frequency" obstruction for the existence of currents satisfying (29). In particular, 
in the Schwarzschild case, in view of the considerations of Section 3.2, one can show 
that there cannot exist a current J'^'™ associated to a vector field X of the form 
(28) satisfying (29) near r = 3M. Cf. classical results due to Ralston [37]. 

Remarkably, it turns out that this trapping obstruction is the only obstruction 
for the validity of (29) away from the horizon, and one can construct a current 
J-''^'"' where X is of the form (28) satisfying 

(31) f K^"-"'[*]>C^fli,i?. f {{r-3Mf\d^f + {dr^f + \^f), 

JB JBn{Ri<r<B.2} 

for all 2A/ < i?i < i?2 < °°- That is to say, apart from the horizon and infinity, the 
control degenerates precisely on the photon sphere r = 3M. 

The vector field X itself vanishes at r = 3M. The construction of the energy cur- 
rent J'^'™ is quite delicate, and the original constructions in fact required spherical 
harmonic decompositions, a requirement that was later overcome. See [8, 14, 16, 7] 
and subsequent papers. 

Applying now (19) in the region described in the previous section to solution tp 
of the wave equation (2), one obtains 

r ( f {{r- Ulfm^ + {dr^f + 1^1^)] dr 

Jti \JT.-rn{R.i<r<R2} I 

(32) <-c..M*L,/i. :L; Ji-iJ-M"-!^ 
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5.3. The red-shift estimate. To properly understand the behaviour of waves on 
the horizon, we liave to take account of the red-shift effect described in Section 3.3. 
In the present section, we define a timelike vector field N whose multiplier current 
captures precisely this effect. This type of current was first discovered in the 
Schwarzschild case in [14], but it turns out, as was shown in [19], that the existence 
of this current is in fact a general property of stationary black hole spacetimes with 
Killing horizons of positive surface gravity. 

In anticipation of later sections, we might as well at the onset give a statement 
applicable to the entire Kerr family, as the Schwarzschild case is no simpler. The- 
orem 7.1 of [19], specialised then to Kerr, yields 

Proposition 5.3.1. Let \a\ < M, ga.M be the Kerr metric and TZ, etc., be as 
before. There exist positive constants b = 6(a,M) and B = B{a,M), parameters 
ri(a,M) > ro(a, M) > r+, and a Lpt-invariant timelike vector field N = N{a,M) on 
TZ such that 

(1) K^[*] > bJ^[^]Nt' for r<ro 

(2) -K^[^'] < B J^[*]iV^, for r>ro 

(3) N = T for r>ri, 

where the currents are defined with respect to QM.a- 

Property 1 above can be thought to embody the exponential damping property 
discussed in Section 3.3. 

Applying (19) with current in the region of Section 5.1, we obtain using 
the above 'red-shift' Proposition 5.3.1 the following estimate for solutions of the 
homogeneous wave equation = on Kerr: 

(33) 




Here r+ + s_ < tq and r+ + s+ > ri . The red-shift effect is responsible for the horizon 
and the space-time term on the left hand-side of (33). It is generally effective only 
in a small neighborhood of the horizon. The term at null infinity appears as a 
consequence oi N = T being time-like near . 

We note that the above estimate degenerates in the extremal case \a\ = M, 
in accordance with the comments at the end of Section 3.3. This is in fact the 
fundamental new difficulty of the extremal problem. See [4]. 

The other manifestation of the red-shift effect is the following commutation state- 
ment 

Proposition 5.3.2. Let g = gM.a and N be as in Proposition 5.3.1. Then, on T-L^ 
with a distinguished generator L, extending L to a translation invariant standard 
null frame Ei, E2, L,Y , then for all k> 

OgiY''-^) = k^Y''^'^^ + ^ c^E['''EpL"''Y"'^^ 

|m|<A:+l,m4<fc 

where Kk > 0. 
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In the Schwarzschild case, Y and T together wiU provide around every point a 
timehke commutator with good properties. This is important for proving pointwise 
bounds. These commutations will also suffice in the \a\ « M case. In the general 
Kerr case, however, besides Y and T, we shall also need to commute with the 
axisymmetric Killing field $ to obtain pointwise bounds, exploiting the fact that 
the span of T and 4> is timelike except at the horizon, where the span of T, $,y 
again contains a timelike direction. 



6. Schwarzschild: Putting it all together 

Wc now proceed to show how the above components can be used to show decay 
type statements, beginning with integrated local energy decay. 

6.1. Integrated local energy decay. The choice of the virial current J"^'™ in 
Section 5.2 is such that the boundary terms 



can be controlled (with the help also of a Hardy inequality) by 



Putting together (32) and the conservation of energy expressed by (22) we thus 
obtain 

ri f {{r-3Mf\dipf + {dr^f + \^pf)]dT 

Jti \Jj:^n{Ri<r<R2} I 

(35) < Cr,,r, f 3l[^]n^. 

From this statement, one can retrieve the usual integrated local decay (away 
from the horizon) with a "loss of a derivative" , namely: 

(36) rV ^ (m' + m')<CR,M2 f (J^V^J + Jp^X. 
or alternatively 

(37) r i ^ (m' + \i^\')<CR,.R2 f (J^V'l + EJn^'^lXr 

It is perhaps worth noting that if Pitp denotes the projection of ip to the space 
of spherical harmonics with number i, the second term on the right hand side of 

(37) can be dominated by the first. Thus, we obtain precisely (35) without the 
degeneration at r = 3M . Fixed spherical harmonics thus do not see the trapping 
obstruction. 

Finally, estimate (35) together with the red-shift estimate (33) yields the non- 
degenerate (cf. the degenerate (25)) energy boundedness statement 

(38) f < C f 
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and 

{{r-3My-\d^\' + idr^f + \^'))dr 

Jti \Jj:rn{2M<r<R2} I 

(39) <Cr, f \dij\\ 

Jti \ JExn{2A/<r<i?2} / 

(40) <Cfl, [ {\d{T^)\^ + \d^f). 

The latter statement is a true integrated local energy decay, without unnecessary 
degeneration at the event horizon. 

6.2. Aside: boundedness without the virial identity. It is a nice exercise for 
the reader to see that (38) can be derived using only the red-shift estimate (33) 
and conservation of energy (22), i.e. without appeal to (35). 

Together with the red-shift commutation Proposition 5.3.2, estimate (38) allows 
one to prove uniform pointwise boundedness for ip and all its derivatives to arbitrary 
order. See [17]. This in particular already retrieves (and strengthens) the celebrated 
Kay-Wald Theorem stated in the introduction. 

6.3. A final ingredient: estimates at null infinity. The final ingredient for 
the statement of decay is a family of estimates derived in the neighborhood of 
future null infinity I^. These estimates have been introduced in [20] and in a 
way serve as a replacement of the older and less flexible approach involving the 
current 3^ with a vectorfield Z constructed as an analog of the conformal vectorficld 
Z = (^2 + 7-2)^ + 2tr-§p of Minkowski space. 

The family of estimates is based on the application of the current with a 
vectorfield V = r^-^ composed from a power of the Schwarzschild r-coordinate and 
the null coordinate vector field -§;;=-§[+ gf? ■ Parameter p is allowed to vary in the 
range < p < 2. The estimates are cleaner when E,- is in fact exactly null in r > i?3 
for some sufficiently large R3. The resulting estimate for solutions ip of the wave 
equation then has the form 

(41) 

[ rP-'\d,,{r^P)\' + [ {prP-'\d,{r4')f + {2-p)rP-'m^) 

<clf y-'\d.iri^)f+ f^' f la^r + H^], 

\Js^-^n{r>R3} Jti JT,rn{R2<r<R3} ) 

and holds with R2 < Rj,. 

6.4. Decay of the energy fiux and pointwise decay. We begin with the case 
of Minkowski space, where already, the p- hierarchy of estimates (41), coupled with 
the statement of integrated local energy decay (30) and the energy conservation 
associated dt , provides a new approach to the problem of robust decay. 

Applying (41) first for p = 2, and using the integrated decay estimate to con- 
trol the second term on the right hand side, one may extract (by the pigeonhole 
principle) a suitable sequence of hypersurfaces T,^ for which the first term on the 
right hand side corresponding to the p = 1 estimate is bounded with n in the role 
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of Ti . Repeating this argument and using also conservation of the usual energy, one 
obtains easily: 

(42) [ < ^ [ r^ldif. 

This leads, after appropriate commutation arguments, to the pointwise decay 

(43) |r^|<-, 

T r 2 

with a constant C depending on the initial data. Sec [20] . 
A further refinement would also give 

(44) H<4. 

r 2 

We note that the estimate (44) matches the strongest rate of decay obtainable by 
the full strength of the usual vector field method. Here, however, we have avoided 
the use of either the conformal multiplier or commuting vector fields generating 
Lorentz rotations and scaling, i.e. multiplier currents or commutators with weights 
in t. This is important for the black hole case, in particular Kerr with \a\ < M 
Moreover, in view of this additional robustness, this new method thus has an even 
a wider range of applications to nonlinear problems. See [48]. 

In the Schwarzschild case, the argument requires in addition commutation with 
T in view of the loss of T on the right hand side of (40). Other than that, it is 
identical to the Minkowski case! Combining (38), (40) and (41) one gets, through 
an analogous iteration process to the above, the decay of non-degenerate energy 

(45) [ |9VI'<4 [ {AdTT^\^ + r''\dTi'\^ + r''\dM^). 

One then obtains (43) by standard commutations. Again, one can then obtain 
further refinements like (44) without introducing i-weiglits. 



7. Kerr: The strategy revisited 

In the Kerr case, we must revisit the Schwarzschild strategy in view of the new 
additional difficulties arising from its more complicated geometry. 

7.1. The difficulties. 



7.1.1. Superradiance. Recall the discussion of the ergoregion S in Section 3.1. For 
p € the energy density J^[^]n''(p) in general fails to be nonnegative. Since if 
a + 0, £ nYiT- t and £ n H'^ + 0, it follows that the first two integrals on the left 
hand side of (22) no longer have positive definite integrands, and (22) no longer 
provides a priori control on 

In particular, the flux of energy to null infinity-the third term on the left hand 
side of (22), which is positive definite-can thus be larger than the initial energy, 
i.e. the term on the right hand side. This is the phenomenon of superradiance, 
first discussed by Zeldovich [49]. Moreover, a priori, as T2 ^ oo, this flux can grow 
without bound. 



'Cf. the loss in the decay rate proportional to a in the original decay result for \a\ « M in [19] . 
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7.1.2. More complicated trapping. We have explained that in the presence of trapped 
null geodesies, considerations of geometric optics necessitate a loss of derivatives 
or a degeneration on the trapped set in the putative integrated local energy decay 
estimate. In Schwarzschild, this corresponds to the photon sphere r = 3M, and 
one could construct a current J'''"'"' where the vector field X vanishes precisely at 
r = 3M such that (31) holds. 

As we have seen in Section 3.2 in the Kerr case for a + the trapped set has 
a more complicated structure, which can be fully explained only in the cotangent 
bundle picture. Moreover, an argument in [2] shows that there does not exist a 
multiplier of type J-^'"^ satisftying an analogue of (31). 

Thus, a classical vector field multiplier to capture trapping is ruled out. 

7.1.3. The coupling. Recall that the estimate (40) is not used alone, but is used in 
conjunction with the conserved energy estimate (22). Thus, even if one can obtain 
a current of the form (28) satisfying a nonnegativity property generalising (31), it 
is not at all clear how one obtains an integrated decay estimate, in view again of 
superradiance, which now manifests in failure to control the boundary terms on the 
right hand side of (40). Thus the two difficulties are in fact coupled. 

7.2. The case \a\ « M. 

7.2.1. Boundedness. The boundedness problem was resolved in [17]. At the heart of 
the argument were the following 2 insights regarding the difficulties of superradiance 
and trapping: 

1. At least formally, using Carter's separation described in the Introduction, we 
may decompose a solution ip into its superradiant part ip^, and its non-superradiant 
part 

(46) V = V'b + V'j, 

i.e. where V'b is supported in the frequency range (5) and in the complement 
of this range. ^ In the case \a\ « M, one easily sees that the superradiant part 
tpii is not trapped, essentially because all future-trapped null geodesies leave the 
ergoregion £. Moreover, one has a more quantitative version of this statement in 
that one can construct a current J'^'"' satisfying an analogue of the non-degenerate 
(29) at least away from the horizon. 

2. One can overcome the difficulty of controlling the boundary terms of a J'''-'"' 
current by the conserved energy, with the help of the red-shift. We first note 
that we may write 




The reader may have already noticed that it is not actually the complete separation which is 
necessary to define this particular type of decomposition, only Fourier analysis in t and if>. See the 
comments at the end of this section. We suppress here the issue of justifying this decomposition 
a priori, since we discuss it in Section 9. 
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where we have used the energy identity for J"^, and the fact that \d^\'^ contain 

the square sum of ah derivatives appearing in the fluxes on their corresponding 
hypersurfaces. The last two terms in the above inequahty appear precisely on the 
left hand side of the red-shift estimate (33). The constant /? can in fact be chosen 
proportional to the value of the specific angular momentum a of the Kerr solution, 
which means that in the slowly rotating case \a\ « M , given an analogue of (29) 
such that (34) is bounded by the right hand side of (47), we may add the analogue 
of (32) and (47) and a small amount of (33), and the estimates close. 

The above is to be applied only to the superradiant part ip^ , as it is only for this 
that we have an analogue of (29). For the non-superradiant part -tpp we simply 
apply the argument of Section 6.2J 

We obtain thus the analogue of (38) for the solution and thus we retrieve the 
Kay-Wald linear stability theorem, but now for Kerr with \a\ « M . 

Let us note finally that the above argument docs not require the complete sepa- 
ration to define (46), but only the Fourier analysis with respect to dt and d^. Thus, 
it can be adapted to the wave equation on more general axisymmetric stationary 
spacetimcs whose metrics are C^'^-close to Schwarzschild, and this is in fact the set- 
ting of [17].^ It is interesting to note that in this very general class, one does not 
have uniqueness of geodesies, and thus, the behaviour of the trapped null orbits 
may be very different from in Schwarzschild. In particular, one does not expect to 
be able to prove quantitative decay results in this class. 

7.2.2. Decay. For the purpose of decay, the estimate (32) is not a tool for obtaining 
(38) but an end in itself. We must thus obtain the analogue of this estimate for 
the whole solution i/j, not just the surperradiant projection ip^. 

As we have seen, for this, in view of [2], one must go beyond classical multiplier 
currents. The property of trapped null geodesies discussed at the end of Section 3.2 
strongly suggests that suitable J'^'*" currents could be constructed if one were 
to microlocalise the construction, to distinguish between geodesies with different 
sets of conserved quantities, and degenerate accordingly at different values of r 
associated to trapped null orbits, all (in the case \a\ « M\) close to 3M. 

Again, Carter's separation^ into modes (4) suggests a particularly natural mi- 
crolocalisation suitable for this problem. In fact, the existence of this separation is 
equivalent to the separability of geodesic flow and thus, it should naturally capture 
properties of geodesic flow. Using this microlocalisation we shall indeed be able 
to construct microfocal versions of the currents J'^'™ ensuring the validity of the 
analogue of (31). We shall see how this is done in Sections 9-11 below. 

The above microlocalisation was first carried out in Chapter 5.3 of [19]. We note 
that, independently, [41] showed a similar estimate in the |a| « M case considered 
here, using a different microlocalisation relying on the general machinery of pseudo- 
differential operators. The integrability of geodesic flow is essential here too. Later, 
a third 'microlocalisation' was achieved in [3] by commutation with a higher order 
differential operator constructed from the Carter tensor. This attractive approach 



In fact, the elementary argument of Section 6.2 was first discovered in the present Kerr 
context [17]. 

^The C^-closeness is what guarantees stabihty of positive-definitivity properties of non- 
degenerate energy currents. 

^Unlike in the case of the boundedness property of Section 7.2.1 above, here, as will be apparent 
we really need the complete separation. 
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allows for a formulation of the whole argument in terms of differential (as opposed 
to pseudodifferential) operators, at the expense of showing results only for higher 
energy norms with strong weights at infinity. 

Given the fact that boundcdness has been proven (cf. Section 7.2.1), one need 
not worry about boundary terms and the difficulty of Section 7.1.3. This yields the 
analogue of the integrated decay statements (39), (40), where the degeneration is 
now on a set containing the physical space projection of the trapped set. Alterna- 
tively, one can apply directly part 2. from Section 7.2.1 to obtain the boundedness 
of the boundary terms, bypassing part 1. As we shall see, however, the insight of 
part 1. of Section 7.2.1 will return with a vengeance in the study of the general 
\a\ < M case! 

7.3. The general case \a\ < M. The nature of geodesic flow as described in Sec- 
tion 3.2 suggests that 'trapping' is qualitatively similar for the whole range of a, 
and thus, microlocalised currents should again be able to capture this difficulty. 
This is indeed the case and is particularly easy to see with respect to the microlo- 
calisation based on Carter's separation, as the latter is so tightly connected with 
the underlying geometry. 

The real difficulty now, however, is trying to apply (47). For in the general case 
\a\ < M, the coupling constant /? is no longer small, and thus the extra terms cannot 
be simply absorbed by adding a small amount of the red-shift estimate. One has 
to identify some non-trivial structure in the coupling of the various components. 

The miracle that saves the day is the following: The insight 1. of Section 7.2.1 
in fact holds for the entire \a\ < M range! That is to say, superradiant frequencies 
are not trapped, in fact, we can quantify this by saying that if the frequency is 
sufficiently large, one can construct a J'^'"'-typc current with K'^ '" sufficiently 
positive to accommodate coupling with (47) for an arbitrary /3. 

This miracle is somewhat unexpected because it is certainly not the case that 
all trapped null geodesies leave the ergoregion £. It would be interesting to under- 
stand better the underlying reasons (if any) behind the phase-space disjointness of 
the trapping and superradiance phenomena. If better understood, in view of the 
programme of [1] , this could have bearing on the uniqueness of Kerr problem. 

The above remarks all concern the situation where (at least some) frequencies are 
large. One must also deal with the case where all frequencies are bounded without 
there being a large frequency parameter to exploit. In the \a\ « M case, these can be 
dealt with by a perturbation argument from Schwarzschild.^*^ Again, however, the 
separation gives a systematic way to deal with these difficulties which is important 
in the large \a\ < M case where perturbation arguments from Schwarzschild are not 
available. See Sections 10 and Section 11.7. 

8. Carter's separation as a geometric microlocalisation 

We have given a complete overview of the proof, and now we begin with the 
technical details. As we have explained in the previous Section, we shall appeal 



This was the approach adopted in [19]. In our alternative approach [18], wc deal with the 
low frequencies as in the present paper (see Section 10), precisely because the method is more 
systematic, can be used in the general |a| < M case, and yields an independent self-contained proof 
of the Schwarzschild result. 
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to Carter's separation of the wave equation, viewing it as a method to frequency- 
locaHze energy estimates in a manner particularly suited to the local and global 
geometry of Kerr. In the present section, we introduce the relevant formalism. 

8.1. Carter's separation. The separation of Hgip = requires taking the Fourier 
transform in t, and then expanding into what are known as oblate spheroidal har- 
monics. As we shall describe in Section 9, we will apply this separation to solutions 
of the inhomogeneous equation 

(48) Ug^ = F, 

where 5* is related to ip by the application of a suitable cutoff, making ^E* then L^{t) 
for each fixed (r, 0,0), where these denote Boycr-Lindquist coordinates on TZ. 
We may thus write '^{t^r,0,(j)) in the form 

r, e, (f)) = —= / e"'"**(w, r, 0, (j))duj 
\/27r J 

and further decompose 

ml 

The collection {S',„£(^, cos0)e*™'^},„£ of the eigenfunctions of the self-adjoint oper- 
ator 

sm oO \ oO I aq)'^ sm 6 
on (sin 9 d9 dcp) forms a complete orthonormal basis on {sin 9 d9 dcj)) . It is 
parametrized by a set of real discrete eigenvalues A^] 

P{0 Smi{tcos9)e'^^ = Xrai{0Sme{^,cos9)e'^t 
which have the property that 

(49) Xmi{0+f>H{\m\ + l). 

For ^ = the oblate spheroidal harmonics S'm£(^,cos0)e""'^ reduce to the standard 
spherical harmonics Ymi with the corresponding eigenvalues A(0)^^ = £{£ + 1). 

The usefulness of the above decomposition in the context of the wave equation 
on Kerr is given by the following result of Carter [10]. If '^{t, r, 9, 0) satisfies 

□gvp = F, 

for g = ga.M a Kerr metric with \a\ < M, with ^'(t,-) in L^{t), then defining the 



coefficients F^^\r) as above, the following holds: 



(50) ={r' + a')AFir\ 

Using the definition (16) of r* and setting 
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we obtain 

(dr*) 



(51) TI^^LT^ + - V^7\r))u - H[:f 



where 



+ 



rai V I {r^ + a^y (r2 + a2)3 (r^ + a2)^ 

In the the Schwarzschild case: 



(52) V^^\r) = {r-2M)i^-^ + ^\ 



^3 ^4 j ' 



■ 2A/ / 2A™<;(3A/ - r) ^ 2M{8M - 3r) 



We note that in this picture the apphcation of ^ to ^' corresponds to the muhi- 
pHcation of ^^l^^ by iuj, while ^ corresponds to the multiphcation by im. 

8.2. Notational conventions. In what follows, let us suppress the dependence of 
u, H and V on auj, m, £ in our notation. 

We will also use the notation ' = -£-^ and introduce a new useful non-negative 
parameter 

which according to (49) obeys 

A > |m|(|TO| + 1). 

A word of warning: For fixed 5j\/,a, we shall often refer to r*-ranges by their 
corresponding r ranges, and functions appearing in most estimates will be written 
in terms of r. Moreover, given an r-paramcter such as R, then R* will denote 
r*{R). It is important to remember at all times that ' always means ^! 

8.3. Boundary conditions. As we shall describe in Section 9, the separation 
process will be applied to v]/ = ^{t)i(; obtained by applying a cut-off ^ with the 
property that ^(r) = 1 for t > 1 and ^(r) = for r < to the solution of the 
homogeneous wave equation ip. As a consequence, the coefficients and 
hence M^^'(r) will satisfy the outgoing boundary conditions: 

f am 

(54) u + i(u} )u = 0, r = r+, 

(55) m' - iciju = 0, r = oo. 
They correspond to the standard outgoing conditions on ^' 

did ad 



[dr* \dt^ 2Mr+d4>)) 



^"'^Before suppressing the dependence, it miglit be useful to remark that in fact, the functions 
u, H, V, etc., depend on both auj and a, whereas the parameters depend only on auj. For the 
former 3, will view reference to the a-dependencc as implicit in the reference to a in auj. Hence, 
our writing Olj instead of in (52). 



THE BLACK HOLE STABILITY PROBLEM FOR LINEAR, SCALAR PERTURBATIONS 25 



We note that the vectorfield L = -g^ + 2j\/r ^ is the distinguished null generator 
of the future horizon "H^ mentioned, for instance, in the earlier discussion of the 
red-shift effect and surface gravity. 

8.4. The separated current templates. To describe the analogue of energy 
currents localized to frequency triplet (cj,m,£), it will be convenient to define the 
following current templates. 

8.4.1. The microlocal 3^ . We first define the microlocal analogue of the conserved 
energy current, i.e. the current: 

^^[u] = cjIto(m'u). 

We compute 

Q'rp = u! lm{Hu). 
We observe that the boundary conditions (54), (55) imply that 

Qrir,) = - ^M', Qt(oo) = io'\u\' 

and that Qx{oo) and -QT{'r+) are both positive in the non super-radiant regime 

a 2M r+ 

8.4.2. The microlocal virial currents J"^'™. Next, we define the microlocal analogue 
of the virial currents J"^'™ where X is in the direction of dr* , and w is a suitable 
function. 

It is most natural to parametrize such choices of currents as follows. For arbitrary 
functions f(r*),h{r*),y{r*), we define the currents 

= / [\uf + (c.2 - V)\u\'] + f'Reiu'u) - i/'>p, 

Q'I = hRe{u'u) - h-i'\uf, 
Qy = y[\uf + {^'-V)\uf]. 
We note the simple relation: Qg = Qi + Qi- We compute 



(QlY = 2/Vf - fV'luf - -f"'\u\^ + 2fRc{u'H) + /'Re(ui7), 



(Q^)' = h [\u'\^ + {V- Lo^)\u\'] - ^h"\uf + hRe{uH), 

(Qir = y' [W? + {^^ - V)\u\^] - yV'\uf + yRc{u'H). 

8.4.3. The microlocal red-shift current . Finally, we define the microlocal ana- 
logue of the red-shift current J^. This takes the form 



Qred 



' ^ 2Mr^' ' ^ ' 2MrJ ' 



where z will be chosen to be a suitable function blowing up as r ^ r+ . 
For convenience we denote 



2Mr+' 

a constant which can be also characterized by the relation 

- V{r+) = \uj - w+p. 
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SO that V{r+) = and V'{r) = V'{r). Wc then compute 

(56) (Qrod)' = z'\u' + - - (zV)' \uf + 2zRe (u'H). 



9. Proof of Theorem 1 
In the rest of this paper, we give detailed account of the proof of Theorem 1. 

9.1. The continuity argument in a. The proof of Theorem 1 proceeds via a 
continuity argument in tlie Kerr parameter a. This is central to the logic of the 
proof because we are relying (recall Section 8.1) on taking the Fourier transform in 
time, and this requires some a priori information. 

We must show thus that the set of a for which the statement of Theorem 1 holds 
(with a constant C = C{a) whose form is given by the proof) is nonempty, and both 
open and closed in the subset topology of {-M,M). Wc note that in principle, the 
function C(a) will diverge as \a\ M. 

The value a = corresponds to Schwarzschild. The validity of the result in the 
Schwarzschild case thus satisfies the requirement on non-emptiness. (In fact, the 
considerations of Section 10 below can be used to provide an independent, self- 
contained proof for the Schwarzschild case. See [18].) 

9.2. Openness. The template for the openness argument is actually the small 
\a\ « M case, treated in detail in [18], which can be thought of as implementing 
the openness argument in the vicinity of a = 0. 

The actual openness argument is an amalgation of techniques from the closedness 
argument and the small a case. For convenience, here we shall treat the openness 
case as if it were exactly the small a case. This will allow us also to review in more 
detail this special case. 

9.2.1. Cutoffs. For openness, wc do not know that ip decays a priori in time, and 
thus, our cutoff ^ must be supported only in < r < for some large Tf, where 
^ = 1 in the range 1 < t < Tf - I say. 

By a bootstrap argument, Tf will finally be taken to infinity. The proximity 
to something known (inherent in proving openness) will be essential here. See 
Section 9.2.3. 

9.2.2. A microlocalised virial current. For convenience, let us actually restrict to the 
|a| « M case in our discussion of this construction, because it is a nice illustration 
of how proximity to a solution where estimates have already been established can 
be used. (The general openness will follow from an amalgamation of the ideas here 
and in the closedness section.) 

We have discussed in Section 7.2.2 that wc need a microlocal version of the 
current J-^-^ and of the positivity statement which schematically we write 



where r+ < Rq < 3M < < oo. The £'-^^'"[*] occurs (recah (19)) from the fact that 
^ = ^-0 and thus satisfies an inhomogeneous wave equation. 
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In the formalism of our separation, this means that for various frequency ranges 
{uj,m,i), by employing various compositions of the virial microlocalised currents 
Qqj Qi I Qfi seek a current Q such that 

Q'[u]> / \uf + x{^ + ujV\ + \uf- / H-{f,h,y)-{uy), 

oo ^0 J —oo 

where the last term is schematic for various combinations (such as fKe{u' H)) which 
occur. The terms Q, u, Xi H, /, h, y should all be understood with the full set of 
indices, i.e. u^'^\ etc., here omitted for brevity. For frequencies for which trapping 
is relevant, x will vanish at a unique value of r associated with the frequency triple. 
One restricts the domain of integration on the first term on the right hand side 
because one expects this virial current not to control things at the horizon and 
infinity. 

We may rewrite (57) above as an estimate 

r r oo oo 

b \u'f + x{X + io')\u'\ + \uf< Q'[u]+ H-{f,h,y)-{u,u') 

•J Rq J — oo J — oo 

and thus 

(58) 6 \u'\^ + x{X + uj^)\u^\ + \u\^ <Q[<^]-Q[-^]+ [°° H-{f,h,y)-{u,u') 

JR* J-oo 

9.2.3. Summing. Restoring now the indices, summing (58) over m and t and inte- 
grating over w, we obtain by Planchcrcl for 7'^ < Rq < Ri < Rc 

(59) + f \dif* f 

Here, x ^lO"^ denotes a smooth cut-off function which vanishes on a small (a-size) 
neighborhood of r = 3M. 

Below we explain the main features of how we obtained (59): 

1. The range of r-intcgration on the left hand side of (72) is restricted to [0,t/]. 
This is due to the fact that our microlocal analysis was applied to the function 
^ = S.{t)'iP with a cut-off function ^ vanishing for r < 0, r > r/. The former 
property of £, led to the microlocal outgoing boundary conditions (54), (55). 

2. The boundary terms on and arise from the corresponding boundary 
terms Q(±oo) in the microlocal virial inequalities. The fact that these terms are 
controlled by the fiuxes is a consequence of the boundary conditions (54), (55). 
For instance, 

Qo(+°°) = + (^^ - V)\uf)i+oo) = 2fuj^\u\^{+oo) = ifuju{-iuju - it')(+oo), 

which by Plancherel can easily be seen to be bounded by the J-^-ffux upon summa- 
tion. 

3. The contribution of the inhomogeneous term H in (58) gives rise to the terms 
integrated over Sq and ^Tf-i on the right hand side of (59). Let us see briefly how 
this occurs. 

Recall that H arises from the term F = ng^t/j + 2X7^^X7^^'; introduced by the 
cut-off £, mutliplied to -0, and that F is thus supported in r e [0, 1] u [r/ - l,r/]. 
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For concreteness, consider a term 



f 



Were this integral in fact restricted to [Ro,Ri], it is easy to see that it can be 
controlled by a little bit of the left hand side of (59) and the last two terms on the 
right hand side. For we apply Cauchy-Schwarz to bound this term by 

2 \ 



/I 



du 



dr* 



Jo JS^r 



and applying Plancherel, one can then absorb the first term into the left hand side, 
while the right hand side, together with a finite time energy estimate and the above 
comment on the support of F, can easily be seen to be bounded by the claimed 
boundary terms. 

On the complement of [i?o,i?i], one needs a slightly different strategy, which 
takes into account that the term H^^\r) arises from the compactly supported 
expression F = □g^(r)5' + 2V^^(t) V^^* with a smooth cut-off function ^(r) = 1 for 
r > 1 and ^(t) = for r < 0. In [18], wc have arranged our mutipliers, i.e., function 
fj^i^r) in this example to be independent of the parameters uj,m,£ outside of 
{Ro,Ri) so that after summation in m,£ and integration in lj we obtain the term 

rn{,i{B^o,Ri)} dr* 

Since F is supported in the interval r € (0,1) we can easily relate the above ex- 
pression to the initial data. More generally, the basic idea for treating compactly 
supported in r inhomogeneous terms in the equation Ogip = F to couple them 
with the based energy estimate for ip. (The precise analog of this for the wave 
equation in Minkowski space is that the integrated local energy decay estimate can 
have two incarnations (with arbitrary a > 0) 

and the second estimate is more advantageous if F has compact support in t.) 

These technical issues are carefully treated in [18] and we proceed without mak- 
ing any further comments. 

9.2.4. Completing the argument. We may now couple (59) to (47), and, applying 
the red-shift estimate, as described in Section 7.2.1, we may obtain 

(60) T'dr f (xiaV-l' + 1^1') < C f 

Jo J'Srn{r<Ri} JSo 

where we are using that we may treat /3 as a small parameter to absorb terms. The 
cutoff time was arbitrary and does not appear in the constant, and thus, we may 
finally replace Tf with oo. 

We have thus reduced openness to showing (58) for all frequency ranges. This 
will be done explicitly in Section 10 for \a\ « M, giving us the opportunity to 
review this simpler case-treated in full detail in [18]-first. As said earlier, the actual 
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openness argument for the general case can be implemented by an amalgamation 
of these multipliers and the types of multipliers we shall use in the closedness case. 
Thus, the real new difficulties of the general \a\ < M case concern closedness. 

9.3. Closedness. For closedness, it suffices to improve an estimate for a solution 
already known to decay, i.e. it suffices to consider a sequence of a for which the 
statement of Theorem 1 holds, and to show that the constant C{ai) of the theorem 
can be uniformly bounded as -> a. 

9.3.1. Cutoffs. In view of the fact that we already know that the solution is in- 
tegrable towards the future, we need only cutoff near the initial data surface Sq. 
That is to say, we apply a cutoff ^ to the solution ?/> of (2) such that ^ = 1 for r > 1 
and ^ = for T < and we then apply the separation of Section %.l to = ^ip. In 
contrast to the openness case, in view of the new cutoff, the inhomogeneous term 

can be thought to arise solely from initial data for on Eq. 

9.3.2. A microlocalised combined current. Whereas in the small \a\ « M case, we 
needed merely to microlocalise the construction of J'''"'"', and the boundary terms 
could be handled using (47) and the smallness of the parameter /3, here, the cou- 
pling between the virial current, the problem of superradiance, and the red-shift is 
nontrivial. To exploit the structure described in Section 7.3, all currents must be 
microlocalised . Essentially, we would like a microlocal version of the inequalities 



L 



(J;f'"' + /3j;r + aJ>^<0, 



where, not only the choice of X but also the coupling constants a, j3 are 'mi- 
crolocalised'. To put it another way, using the linearity of the currents, we seek a 
microlocalised version of the combined currents: 



J. 



X+PN+aT,w 



In our formalism (see Section 8.4), what we are seeking, for each frequency triple 
(cj,m,^), a combination Q of Qt, Q^cd ^^"^ ^^^^ virial currents Qq, Qi, Q2 such that 
(57) holds and moreover 

Q[oo] - Q[-oo] > 
so we may drop the boundary terms and write 

(61) b f''"\uf + x{>^ + ^^)\u? + \u?< H-{f,h,y)-{u,u') . 

9.3.3. Summing and completing the argument. If we succeed in obtaining (61) for 
all frequency triples, then as in Section 9.2.3, we may sum, where now we obtain 
simply: 

(62) 6 Pdr r {x\d^\U\dr^^\' + \4>\')< f \di;\\ 

The difficulty of the integral does not occur now in view of the new definition 
of the cutoffs. The rest of the proof follows as in Section 9.2.4, completing the 
argument for closedness. 

We have thus reduced closedness to showing (61) for all frequency ranges. We 
shall show this in Section 11, completing thus the proof of Theorem 1. 
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10. The SMALL a case 

This section and the next will contain the microlocalised estimates appealed to 
in the proof of Theorem 1 as given in Section 9. We begin by partitioning the 
frequency space into relevant frequency ranges. 

10.1. The frequency ranges. Let uji, Ai be (potentially large) parameters to be 
determined, and A2 be a (potentially small) parameter to be determined. We define 
the frequency ranges T^, T-^, JF^, by 

• = {(a;,m,£) : |a;| < wi, A,„f(aw) < Ai} 

• = {(w,m,£) : |a;| < ui, \,nt{auj) > Ai} 

• = {{uj,m,£) : |a;| >a;i, Xrni{aui) > A2W^} 

• = {{uj,m,e) : \uj\ > wi, X„ie{auj) < A2W^}. 

Our constructions will have a separate philosophy for each of the above ranges. 

10.2. The JTjj range (bounded frequencies). This is a compact frequency range 
and, in view of the fact that we have already constructed multipliers in the Schwarz- 
schild case (see [14, 16]), by stability considerations, their positivity properties carry 
over for \a\ « M. See our proof in [19]. 

In [18], we chose to give a different argument, which, on one hand, allowed us 
to treat the general axi-symmetric case (jaj < M,m = 0), and furthermore, is more 
useful in applications to the openness regime of the general case. 

Let us define 

AMramu) - a^m? 



A 



(r2 + a2)4 



(r2 + a2)2 

((A„^ + w2a2)(r2 + a^f + {2Mr^ + a^r^ + a" - AMra^)) . 



(63) 

We have 

y+ > 

for r > Moreover, according to our conventions, for all jaj < « M, 

B{A/r^)Ar-^ + Br'^ >V+> {A/r^)bAr-^ + b{A/r^)r-^. 

We may now fix an arbitrary r+ + s_ > Tc > 2M, where s_ is as in (33) and can be 
chosen only to depend on as, and it follows that for sufficiently small depending 
in particular on this choice, we have for jaj < and (a;,m,£) e J^^ the following 
inequality in the region r > rc'. 

(64) B{A/r'^)Ar-^ + Br'^ >V> {A/r^)bAr-^ + b{A/r^)r-^, 
whereas, for all r > r+, we have 

(65) \V'{r*)\ < B(A/r^)Ar-^ + B{A/r^)r-\ 

Finally, again choosing as sufficiently small, it follows that for all r < rc 

(66) y'(r*) >6(A/r2)A, 
whereas there exists a constant Rc such that for r > Rc, 

(67) -V'{r*) > {Alr^)bAr-^ + b{A/r^)r-^ 
in this frequency range. 
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We will now split the frequency range into two subcases, considering each 
separately. 

10.2.1. The subrange \u)\ < ljt, (the near- stationary subcase). The motivation for the 
current to be constructed here is that in the Schwarzschild or m = case, applying 
Qi with h = 1 immediately excludes nontrivial stationary solutions w = 0. 

We will fix an cl>3 > which will be constrained in this subsection to be small. 
Because uj^ is not exactly the naive current h = 1 must be modified. 

We begin with a current which will be defined with an h = h{r*). 

Note first that given arbitrary oq < M and q > 0, p > 0, i?c > such that e^^ Rc 
is sufficiently large and p sufficiently small, for each \a\ < Os, we can define a function 
h{r*), such that the following hold: For r < rc, 

Q<h<R-\ \h"{r*)\<ql\r*\\ 

with h{r*) moreover of compact support when restricted to -oo < r* < r*, whereas 
for rc<r< e'P ^ R^ 

h>R-^-A/r^, h"{r*)<0, 
whereas for e"^ Rc < r < Rc, 

\h'{r*)\ < 4Rfp/r, \h"{r*)\ < 4R-%/r^, 
whereas, finally, for r > Rc, 

h = 0. 

This h will be useful in view of the positivity of V+ . We note finally that h can be 
chosen so that h restricted to r > rc is independent of the parameter q. 

Let us also consider a current Q2 defined with a y = y{r*). Like h, the function y 
will be independent of uj, m, and £ in the allowed range. Given the parameters Rc, 
p, with e^P i?3 sufficiently large and p sufficiently small, the function will satisfy 
the following properties: We set y(r*) = 0, and for r < rc, 

y'{r*) = h, 

noting that y is bounded below in r < rc by a potentially large negative constant 
depending on q, in view of the fact that h is identically for sufficiently low r* . 
For rc < r < e^^ Rc, we require 

yV)>0, {yVy{r*)<lR:'Nr', 

whereas for e"^ Rc < r < Rc 

y'{r*)>bRf, -{yVy{r*)>bRf/r', 
whereas finally, for r > Rc, 

y=l. 

In general, for sufficiently large e^^ Rc and small p, we can indeed construct such 
a ?/ in view essentially of (64), (65) and (67). Let us add that y restricted to r > rc 
can be chosen independently of the choice of h. 

Consider now the current Q = Qi + Q^- In r < rc, recall the one-sided bound 

-yV > 

which follows from (66). It follows that in r < rc, 

Q'{r*) > -g|r*r>p + hRe{uH) + 2yRe{u'H). 
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For < r < e P i?c, choosing p appropriately, choosing ^3, sufficiently small 

Q'ir") > R-^{A/r^)\u'\^ + bV^uf + hRc{uH) + 2yRc{u'H). 
For e"^ Rc < r < Rc, we have for p, uj^ suitably small 

Q'{r*) > bRf{A/r^)\u'f + hRe{uH) + 2yRe{u'H). 



Finally, for r > Rc wc have, 



Q'{r*) = {Qyy{n>2yRe{u'H). 



We obtain finally. 



fe(Ai) r^"(A/r2)|uf + (A/r2)r-3(l + A + oj^)\uf dr* 




(68) </ " q\r*\-^\ufdr* + / {2yRe{u' H) + hRc{Hu)) dr* 




+ Q|(oo)-(Qj+Q|)(-oo). 



The above is the prototype of the type of inequality we derive for all frequencies. 
Cf. (58). Note, however, the extra term proportional to q. We shall deal with this 
in Section 10.6. 

10.2.2. The subrange \uj\ > W3 (the non- stationary subcase). The construction of 
this section will yield a positive current for u)^ small and uji, Ai large depending on 
the smallness of as- 

Let us note that in fact, the construction of this section is quite general and 
depends only on the (almost) positivity properties of the potential V together with 
its asymptotic properties, and not, in particular, on the sign of V' , which will be 
crucial in Section 10.4. We note that similar constructions have a long tradition 
in spectral theory and arc typically used to prove continuity of the spectrum away 
from the uj = 0. 

We consider a Q2"Current. Recall that we have 





< 



00 00 

/ y'\uf + 2l {y'v:/y')\uf 



J 00 ^ — 00 



+ 2/14|M|2(oo)-yKH2(-oo) 




We now define 



(71) 



y = e 



2aj3^ /~ V+dr* 



Note that in this case. 



we have that 



y' = 2uj^^V+e 



2u]3^ V+dr* 
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We thus have for this choice of y that 

/oo 1 oo 1 oo 

{yv.)'\u\^ < i / y'\uf + ^ f u^lvW- 
oo Z J — oo Z J —oo 

The potential Vqcw obeys the estimates 

Therefore, the term ^ujly'\u\'^ dominates both -y'Vncw\u\'^ and -yV^'g^lup in (70), 
provided that « W3. Integrating, we finaUy obtain the inequahty: 

6(c.3,c.i, Ai) [''\A/r')\uf + (A/r2)(l + {A + l,^)\u\' 

/OO _ 
2yRc{u'H) + Qlioo) - Qli-oo). 
00 

10.3. The range (angular dominated frequencies). The considerations 
of this section will constrain Ai to be suitably large, depending on the choice of 
wi. This range is in fact covered by the arguments in Section 11 for the general 
case \a\ < M, and, of course, in [18]. Briefly, a simple construction in this regime 
involves the current Q = + Qi with a positive and monotonically increasing / 
and a function h of compact support peaking in the region where V dominates . 
The final resulting estimate is: 

b f^^ (|uf + (A/r2)r-^(l + (A + u^)\uf) 

Jr* 

< f°° 2fRc{u'H) + {f' + h)Rc{Hu) 

J — 00 

+ Qi{^)-Qi{-°°)- 

10.4. The JF^ range (trapped frequencies). This is the frequency range of trap- 
ping. It is once again subsumed by the trapped case Q'^ which is treated in detail 
in Section 11. 

Thus, we take the opportunity here to merely highlight the simpler construction 
in [18] which is special to the |a| « M case. This is a good warm-up for following 
the detailed but more technical construction of Q'^ in Section 11.5. 

The main observation is that V' has a unique non- degenerate zero r)^^ for 

(uj,m,i) 6 J^^, depending smoothly on w. This can be thought of as a quantised 
version of the property of trapped null orbits described at the end of Section 3.2. 
We consider the current Q^: 

{Qir = 2/Vf - fV'\u\' - + 2fRe{u'H) + f'Rc{uH). 

The function / is chosen in such a way that 

(1) /V)>o, 

(2) / < for r < rl^f and / > for r > r^^f, 

(3) /"' < in the neighborhood of rl^f, 

(4) -fV > dominates outside of a small neighborhood r^^'^K 
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We obtain finally 

6^^'^^((A/r2)r-2|u'|2 + (A/r2)r-2|u|2 + A/^V5(r-rJ,7))2(A + L^2)|^|2j 
< f°° 2,fRe{u'H) + f'Re{Hu) + Q^(oo) - g^(-oo). 

^ — oo 

10.5. The range (time-dominated frequencies). Again, this range is sub- 
sumed by the regime in the general \a\ < M treated in detail in Section 11, so 
here we can be brief. 

For small enough A2 and large enough ui, we have 

2 

for frequencies in . We apply the current Q% with a function y chosen so that 
y>0.,y'> B{A/r'^)r^'^ in (r+, and j/ = 1 for y > i?c + 1. Since, for (w, m,£) e 
we have 

\V'\<B{A/r^)r-\{X^e + aluj^) + l). 



We obtain 



b f''^A/r\~^{\uf + {u;^+A)\u\^) 

Jr* 

< p2yRe(u'5)+g|(oo)-Q|(-oo). 

^ — 00 

See Section 11 (or alternatively [18]) for more details. 

10.6. A note on summation and the red-shift. We have obtained inequalities 
(58) for all frequency triples, except for an extra term 

['\\r*\-^\u\''dr* 
J —00 

coming from (68) corresponding to the < UJ3 sub-case of the J-^ frequency range. 
After summing, we thus in fact obtain: 

(72) 

J^rn{re[Ra,Ri]} JSq Jq JT,rn{r<ra} 

We now add the red-shift estimate (33) and the 3^ identity (47) to eliminate 
the small g-term and the contributions. We obtain again (cf (60)) 

b f^dr f (x|9^|2 + H')< f 

Jo JSrn{re[Ro,Ri]} J^o 

This is the desired integrated local energy decay estimate, which may now be cou- 
pled with the red-shift estimate and the estimates at null infinity (41) of section 
5 to obtain the statements of decay for the non-degenerate energy and pointwise 
decay for the solution of Theorem 1 . 
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11. The GENERAL \a\ < M case 



In Section 9.3.3, we have reduced the problem to obtaining the inequality (61) 
for all frequency ranges. In this section, we shall obtain precisely these inequalities, 
completing thus the proof. 

We begin with a closer examination of some properties of the potential V = 
V^^\r). Without loss of generality we may assume that the specific angular 
momentum a is positive < a < M, as the problem admits a symmetry with respect 
to simultaneous reversal of the signs of a and the azimuthal frequency m. 



11.1. Properties of the potential V . We decompose the potential V =Vq 
with 



A(3r2 -4Mr + a2) i/^^r"^ 



(r2 + a2)3 (r2+a2)4- 
We note that Vi does not contain any frequency parameters and that 



(r2 +a2)4 

Lemma 11.1.1. For any values of the parameters A > 0, m and u the potential 
function Vq on the interval (r+,oo) has at most one maximum and one minimum, 
achieved at the respective points r^-^^^ and r^^^^ with the property that rJJ^jn < f maa; • 
Moreover, for all sufficiently large A the value rjjj^^ is hounded uniformly from above 
provided that either muj > or a^uj^ < CA for some constant C , in which case the 
hound for r'^^^ may depend on C. 

Proof. We compute 



— Vo = 4maMuj 
dr 



4r^ 



+ a 



2)2 (r2 + a2) 



2\3 



Ara'^m'^ 



A 



(r2 + a2)3 (r2 + a2^2 



2(r - M) 



4rA 



1 


(r2 


+ a2)3 




A 


(r2 


+ a2)3 




1 


(r2 


+ a2)3 



Thus, wc have 
d 
dr 



[AmaMu{-3r^ + a^) + Ara'^m'^] 



{2r^ + 2aV - 2Mr'^ - 2Ma'^ - ir^ + 8Mr^ - ira^) 



[4maMa;(-3r^ + a^) + Ara'^m^ - 2A{r^ + a^r - SMr'^ + Ma^)] , 

6Mr + a^) 



ar 



-24Mamu!r + 4a m 



-6A 



2A(3r- 

„2 



where we set 



2Mr + 4:Mra + 



amuj 



-a — — 

3 A 



A 

We now examine the points where ^ \_{r^ + 1^)''^^]] vanishes. Such points are 
solutions of the equation 



2A/r(l-2cr) + —(1 
3 



2m" 



) = 0, 
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let us denote them by: 



= M(l - 2a) - loY - ^(1 - 

We recall the horizon is determined by the value r+ : r^-2Mr++a^ = 0. In particular, 
r+ > M for all values of a < M. 

If cr > then r2 < M and the only point where ^ [(r^ + a^)^^Vi3] = on the 
interval (r+, oo) is 



n = M(l - 2a) + ^M2(l - 2a)2 - ^(1 - ^). 

In that case we note that the function (r^ + a^)'^^Vb^-oo asr^oo, as long as 
A> 0. Therefore, j^V^ can vanish at most at two points r\^^ and r^raax^ where 



j-1 n rl ^ 'T'O 

cxina J77,2J2 ^ max' 

If ri is not real then [(r^ + a^f-^Vo] remains negative for all r > r+, which 
means that can only vanish at one point, where Vb has a maximum. 
If (T < then we must examine 



r2 = M(l-2a)-\/il/2(l - 2a)2 - ^(l - ^) = M(l-2a) 



Since a < M and cr < 0, we have 



\| 3M2(l-2a)2 



Wc note that for < x < ^ , 



which implies that 



3M2(l-2cr)2 3 



2a; 



Vl-a;> 1- 

3 



2M(l-2a)a2(l-2^) 2a2(l-2^) 

To < = < -M. 

9Af2(l-2cr)2 9M(l-2cr) 9 

This now implies that ri is the only zero 

of [(^^ + «^)^i^^o] on the interval 

[r+,oo) and the previous argument applies. 

The last statement of the lemma easily follows from observing that for all suf- 
ficiently large A the leading order behavior of [(r^ + a^)'^^Vb] as r ^ oo is deter- 
mined by the expression 

{6AM - 12Mamuj)r^ - 2Ar^ 

□ 

The next statement effectively establishes that even if rmin exists, it can only be 
'trapped' for the value u = oj+. 
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Lemma 11.1.2. For all values of the parameters A,uj,m, we have that on the 
horizon 

uj^ > y(r+) 

with equality achieved only for lo = In particular, in the notations of the previous 
lemma, this implies that 

OJ^ > Vb(r„„J. 

Proof. Wc compute 

- y = ^2 — — 

{r%_ + a^Y 

1 r < , r9 9 9 91 (2Mr+uj - am)^ 

[AM^rXuj - AMr+amuj + a^m^\- ^ ' 



AA'Prl L + + J UPrl 

The equality 

oj^ = V{r^) 

occurs precisely at the threshold of the super-radiance condition 

am 



□ 



2Mr+ 

We now argue that in the super-radiant regime Vb has only a maximum, that is the 
point rmin is absent. 

Lemma 11.1.3. Let 

am^ 

mui < . 

2Mr+ 

Then 

^V{r^)>^Vo{r^)>cA>0. 
ar ar 

Since we have already shown that r^^^^, if exists, has the property that < r"^^^ 
this implies that r'^^^ does not exist and the potential Vq has its unique critical point 
at rl 



max 



Remark 11.1. The conclusion of the lemma can he easily extended to the range 

2 9 

am . 

amuj < + aA 

2Mr+ 

for a sufficiently small constant a. The smallness of a depends on the smallness of 
the quantity (M - a) . 

Proof. Wc compute 

d AmaMu r o 2] Ar+a'^m'^ 2{r+-M)K 

^^0^^-^ = ^rl + a^Y ^ ^ " ^ ^ {rl + a^f ^ {rl + a^Y 

= , 9 \AmaMuj\-irl + a^l + Ar+a^m^ + 2(rl + a'^)(r+ - M)A] 

(r+ + a-')-^ 

For the values of the parameters mw < the conclusion of the lemma is obvious. 
Otherwise, using the condition 

am^ 

muj < 



2Mr+ 
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we obtain 



{ri + ay-V,{r,)> 
ar 



■{-irl + a^) + Ar+a^m^ + 2{rl + a^){r+ - M)A 



{-rl + a^) + 2{rl + a^){r+ - M)A 



= 2(r+ - M) [A{rl + a^) - 2a^m^] 
= 4(r+ - M) [AMr+ - aV^] . 

The inequalities A > and r+ > M > a imply that ^Vo(r+) > 0. We finish the 
proof by observing that 

dr (ri + a^)^ 

□ 

Our final result of this section gives a mathematical embodiment to the miracle 
which is key to the validity of Theorem 1, namely, the fact that superradiant 
frequencies are not 'trapped'. 

Lemma 11.1.4. For the values < a < M and 



< muj < 



we have 



Proof. Let ro > r+ be such that 



am' 
2Mro ■ 



We compute 



- Vo{ro) = uj 



2 AMroamu) - a^rr? + A A 



1 



1 



{r\ + a2)2 

[(ro + a^)^a;^ - 4AfroaTOw + a^m^ - A A] 



\M\i^r1ij? - AMroamuj + (I'rr? + 



Lj2(r2 - 2A/ro + a^){rl + 2A/ro + a^) - AA 



((r2 + a2)2-4Mro2)-AA] 



A 



9 9 

am 
4M2 



(r^ + a2)2 
(1 + + — )-A 



We now recall that a< M <rQ and that A > m{m + 1) to conclude that 

- K)(ro) < 0. 



□ 
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11.2. The frequency ranges. Let wi, Ai be (potentially large) parameters to be 
determined, and A2 be a (potentially small) parameter to be determined. We define 
the frequency ranges G^, G[^, Gp 5" by 

. ett={(c.,m,0 : c.2 + A>Ai,TOL.6[0,^]}, 

. G^={{co,m,e) : |c.|>c.i, A<A2w2,mc.j^[0,^]}, 

. G^ = {{i^,m,e) : \u\>uj,, A2A>c.2,mc.^[0,^]}, 

. G^^ = {{uJ,m,e) : |c.|>c.i, A2A<c.2<A2iA,mu;^[0,^]}, 
• G\, = {{uj,m,e) : \uj\<uji, A< Ai}. 

11.3. The t/tl range. This a large frequency super-radiant regime. The lower 
bound A > 771(771 +1) together with the superradiant condition implies that in this 
regime 

\ 2A/r+ ) 

and therefore A is necessarily large. In the superradiant range the conclusions of 
Lemma 11.1.3 and Lemma 11.1.4 apply. In particular, the potential Vq has its only 
critical point-a maximum-at r^^^^ and ui'^ < Vo{r^-,^^^). The latter in fact means 
that 

Vo{r) - > cA, V r € [r,nax - S, r„iax + S] 
for some positive constant c and 6 > 0. Note that S can be chosen to be universal 
and independent of the size of (large) A. We further observe that r^jjj. is uniformly 
bounded and is away from the value of r+ : 

c< {rniax -r+) <C 

with the constants c, C independent of A. We also have that the full potential 
V = Vq + Vi has the property that 

±V{r,)>^Vo{r,)>cA. 
ar dr 

Wc now argue that the properties of V mimic precisely those of Vq with the value 
"^max replaced by rraax- 

Vmax fmaxl - . 

In the regime under consideration, the second derivative 

may only change its sign once at a point ri < rJJ^^^ . If it happens then 

^Vo{r)>cA, 
dr 

perhaps with a slightly different value of c for all r e [r+,ri] and. since 



<CA, 



we have that 

ir°nax-ri) > c. 

On the interval [ri, 00), we have 

d { , 9^qCi 

dr 



■{ir' + ayfvoir))<0, 
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and we can choose a value r'^ e (j"!, ^'maa;) such that 



|:Fo(r)>^A, Vr6K,r'i] 



and 



— ({r^ + a^f — V^ir)] < -cAr^ Vr e [r^,oo). 
dr \ dr I 

On the other hand, the potential Vi obeys 

\Vi\<Cr-\ \±V,{r)\<Cr-\ \± i^r^ + a'f ^V,{r)\\< Cr. 

dr dr \ dr I 

Therefore, V cannot not have any critical points on [r+,r5^] and has a unique 
maximum rmax s The conclusion about the distance between r^ax and 

''maa; HOW equally straightforward. 

We may now deduce the existence of an interval {r„iax ~ <5, r^ax + with the 
property that 

V{r)-uj^>cA, ^re{r„,ax-S, ' max 

+ S) 

and 

(73) -(r - r^ax)^V{r) > c(J)A ^'' " ^7°"^' , Vr6[r„oo). 

dr 

We now proceed to the construction of a suitable current for the regime The 
current will be of the form: 

Q = Q{, + Ql + rQ:.,^-BQT. 
To elucidate the exposition we describe the procedure in stages. 

1. We first apply current with / such that / = -latr = r+,/ = Oatr = rmax> 
/ = 1 at r = oo and f'{r) > for all r, giving 

[2/'|uf - fV'luf - \nuf] = [\uf + - + [\uf + 

[2fRc{u'H) + f'Rc{uH)]. 
We impose additional conditions on the function / to guarantee that the expression 

-fV'-^f"'>0 

for all values of r e [r+, oo). Since / vanishes at r = r^ax and V' obeys the property 
(73), we can easily arrange 

It remains to choose / with the property that f"'{r) < in a small neighborhood 
ofr„a, and \f"{r)\<CAr~\ 

The left hand side of (74) is non- negative yet degenerate at r = rmax- The 
condition V{rmax)-^^ ^ cA indicates this regime is non-trapped and the degeneracy 
may be removed with the help of the current Q\. The problem however is a lack 
of control of the boundary terms on the right hand side due to the super-radiant 
condition. The latter will be handled with the help of the red-shift current Q^d- 
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2. We now add a large multiple of the Qi current with a non- negative function 
h supported in [rmax - S,rmax + S] (that is in the region where V{r) > w^), equal 
to one on [rmax - 5/2, 

' max 

+ S/2]: 

1 

2' 



(75) 



(2/' + Ah)\uf + {Ah{V - cu^) - fV) \uf - ^if" + Ah")\u\^ 

Jf~- 00 
[2fRe{u'H) + (/' + Ah) Rc{uH)] . 



The left hand side of (75) is now positive and non-degenerate. Moreover, its inte- 
grand has the property that it is 

S S 

> A(|wf + A|w|^), V r € [r^ax - -, rmax + -]• 

The size of constant A is limited only by the lower bound on the (large) frequency 
parameter A. 

3. We now choose another large constant F with the property that T « A and 
add a F multiple of the red-shift current Q^^^ with z = A-^ and ^ = 1 on [r+,rmaa;] 
and vanishing for r > r^ax + S/2: 



■f 



A 



CV, . 



(^) 



\u' + i{uj - uj+)u\'' - C'\u\^ 



[(2/' + Ah)\uf + {AhiV - u;') - fV) \u\^ - i(/"' + Ah") 

= [\u'\' +((..- ,.,)2 _ rA) + [\u'\^ + 



2(/ - A-i) Rc{u'H) + (/' + Ah) Rc{uH) 



The first term on the left hand side above is positive. To estimate the second term 
we note that we may assume that on [rmax T'''max + S/2], the potential V has a lower 
bound cA. Therefore, its integrand can be estimated by 

< FC [|uf + (cj2 + m^)\uf] « A [|uf + A|up] , 

where the replacement of + m? by A is justified by the superradiant condition. 
Therefore the left hand side is positive definite. On the other hand, using the 
boundary conditions (54), the boundary term at r = r+ on the right hand side is 

{2{u-co,)'-rA)\uf<-^A\uf 

and thus now has a "good" sign and a large multiple constant. The other boundary 
term at r = 00 still has the "wrong" sign but it is merely 

2w>|2 

and therefore can be eliminated by a subtraction of a ^ multiple of the Qt current 

Im(u7J), 

»•+ 

which will create a bounded negative term at r = r+. The latter in turn is easily 
dominated by ^A|up, provided that B « F. We can summarize the results obtained 
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ill this regime in the following inequality: 



A 



r 



f 



(c.2 + A)[H2|__ + l 



[°° -2(/- A^)Re(u'H) - {f + Ah)Rc{uH) + Bulm{uH) 



To end this section we observe that aceording to Remark 11.1 the frequency range 
G^ can be slightly extended to 



< muj < + aA 

2Mr. 



for some small constant A. This observation, in particular, will allow us to restrict 
the range to 



1^ 

2Mr+ 



. ^ am , 
mto f [0, + aAJ. 



11.4. The f/jj range. This is a large frequency non-superradiant regime in which 
Lo dominates all the other frequencies, that is A and consequently to. Since this a 
non-superradiant range, the boundary terms at r = r+ and r = oo of all currents are 
dominated by the (positive definite) boundary terms of the Qt current. 

In this regime, we shall apply a current 

Q = Qi- Qi - bQt 

with a monotonically increasing positive / satisfying 1 < / < 2, and the function 
h = f. We obtain 

Jf oo 
+ (/V-^)-/nH'] 

- [\uf - ((^ -u:,r + Bciu; - u,)) \uW^^^ - [\uf + (1 - B)uM\.^ 

Jf- oo 
[2fRe{u'H) + Bujlm{uH)] . 

Given that > CA with a sufficiently large constant C and 

\V\ < CiA + ciij'^ , ll^'l < ^(CiA + cia;2) 

for some universal constants (large Ci and small ci, which can be chosen to be 
arbitrarily small depending on th size of Ci and a lower bound on > wi) ), we 
may choose C so that the integrand on the left hand side is positive and dominates 

c,f'[\uf + {u;' + A)\uf]. 

The boundary terms are positive due to the boundary conditions u' = -iLd{uj - uj+) 
at r = r+ and u' = iw at r = oo, the non-supcrradiance condition, and a choice of a 
sufficiently large constant B. Therefore, with an appropriate choice of / 

^ - [\uf + (C.2 + A)\uf] + (C.2 + A) [\u\\^r^^^ + \u\\-.rj) 

/- OO 

< / [-2fRc{u'H) + Bujlm{uH)] . 
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11.5. The range. This is a non-superradiant regime in which o;^ is comparable 
to A. This range corresponds to the trapping regime. For these frequencies potential 
Vq may have at most two critical points < rj^^^. 

According to the Remark 11.1 and the discussion at the end of section 11.3 we 
may in fact assume that mo; ^ [0, ^""^ +aA] for some small constant a. Under this 
condition and since lo ~ A, 

io^ -V{r+) =uj^ -Vo{r+)>cA. 
We define tq to be the largest value with the property that for all r e [r+ , rg ) 

Voir) < Vo{r^) + ^A. 

If ro is finite then, since < r*^^^ , it is clear that if r^^^ exists then 

r° < rn < r° 

rrnn ^ ' v — ' max 

Moreover, in that case, since Vb(r) ^ as r ^ oo, the potential Vb does have a 
maximum point r^^^^ and the value r^^^ is uniformly (in A) bounded from above 
with a bound dependent on A2. On the other hand, since |^Vb(?')| < CAr^^, the 
value To, and thus r'^^^, is uniformly bounded away from r+. 
We now consider the value ri , where the function 

d 

tr^ + n'-)'' — Vntr) = -(iA r^ - '/ M r + 4 III rn + 

3 3 A 



dr 



(r' + a^^Voir) 
dr 



r^ -2Mr + AMra + 



changes its sign from positive to negative. If ri belongs to the interval [r+, 00), then 
ri < rj^oj,. Moreover, since in this regime the parameter a = amui/A is bounded, 
the value of ri is uniformly bounded from above. We consider two cases based on 
the value of Vo{ri). 

If Vo(ri) < Vo{r+) + ^A, then, in view of the fact that Vb has a unique maximum 
at r^^^, we have that 

3c 

Voir) < Vo{r+) + —A, Vr e [r+,max(ro,ri)]. 

Moreover, we can find a small constant 6 > such that the above property holds 
on a slightly larger interval 

3c 

Voir) < Vo{r+) + —A, Vr 6 [r+,max(ro,ri) + (5], 



and 

±\ir',ar-^Vo{r) 
dr I dr 



-c((5)A, Vr e [max(ro, ri) + 5, 00). 



These conclusions arc consequences of the simple bound |^Vb(r)| < CAr and 
the fact that A"^^ [(r^ + a^)'^^Vb(r)] is a quadratic polynomial with bounded 
coefficients vanishing at the unique point ri on the interval [r+,00). 

If, on the other hand, Vb(ri) > Vb(r+) + ^A, then we have that (ri - ro) is 
uniformly bounded from below by a small positive constant Ci, and since ^Vb(ro) > 
0, we can find a value r'^ e [ro, ri] such that 

Vo{r)<Vo{r^) + jA, Vr€K,r[,] 

and 

— Fo(r-)>ciA, Vr6[ro,ri]. 
dr 
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Moreover, the last property can be easily extended to a slightly larger interval 



-Voir) > ciA, 
ar 



Vr 6 [ro,ri +6], 



so that, in addition 



d 
dr 



(r' + ay-Voir) 
dr 



<-c{S)A, Vr 6 [ri + (5, oo). 



In both cases, we may therefore claim the existence of the value ^3, uniformly 
bounded away from r+ , such that 



Fo(r)<K)(r,) + ^A<c.2-£A, 
and, such that, for any re [ra, 00), either 



Vr € [r+,r3] 



— Voir)>ciA 
dr 



d_ 

dr 



{r'+ay^Vo{r) 
dr 



< -c{5)K. 



We note that rs is either uniformly bounded from above or rs = 00. 
Now adding a bounded potential Vi we obtain that for r e [r+,r3] 



' I*; 

while on the interval [ra, 00), if non-empty, potential V has a unique non-degenerate 
maximum at r,nax, A^^-close to r^^^, and 

^V{rmax) < -ciA, ^V{r^) > ciA. 
dr^ dr 

In this regime we use a combination 

Q = Qi-Ql- bQt 

of the currents Qo,Q2 and Qt- The current Qq is applied with a function / which 
vanishes on the horizon, monotonically increases on [r3,oo), changes its sign from 
negative to positive at r = r„iax and approaches 1 as r ^ 00, and satisfies the 
additional condition that 

/"'(r)<0, Vr€[r3,oo) 

The second current will be Q^, restricted to the interval [r+,r3], with a monotoni- 
cally decreasing function y. Finally, we subtract a large multiple of the current Qt- 
We obtain for Q = Q{, - - BQt: 

1 



2f'\uf-{fV' + -r)\u\' 



- £ [y' {\uf + (a;^ - VM) - yV'\u\^] + £ 

+ [-f\uf + {B- + [Bujiuj - CO,) - 2y{L0 - LO.f] \u\%-_r, 

= - f [2fRc{u'H)- f'Re{uH)-Bujlm{uH)]+ f ^ [yRe{uH)] 

By the described properties of the potential V, the expression -{fV + \f"') is 
positive on the interval [r3, 00). On the interval [r+,r3], we choose a function y so 
that 

-y'{u:^-V)+yV'-UV'+]-r)>()- 
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Since for these values of r 

,2 



\V'\<CK^, |/| + |/"'|<c4 



for sonic positive constants ci,C, it suffices to to fulfill the inequality 

-±y>yC+C, 

ar 

provided that the constructed function y is positive and monotonically decreasing. 
The function 



Jr 



satisfies all the above criteria. Finally, a choice of a large constant positive constant 
B together with the non super-radiant condition muj ^ [0,7Tiaj+] and the boundary 
condition u' = iuiu at r = oo ensure that both boundary terms at r = r+ and r = oo 
are positive. 

We summarise the results of this regime with the following inequality, which 
holds under appropriate assumptions on /: 



,'|2 , /, ,2 , A \ ( 



^ "^tnax ) I 1 2 



■1) 

[-2fRc{u'H) + f'Re{uH) + Bujlm{uH)] + / [yRc{u'H)] 

r+ Jr+ 



11.6. The range. This is a non-superradiant regime in which the frequency 
is much smaller than A. In particular, by Remark 11.1, the potential Vq has only 
one critical point-a maximum-at the value r = r^ax- Repeating the arguments 
which we applied to the regime we also conclude that the potential V has a 
unique non-degenerate critical maximal point at rmax , which is uniformly bounded 
away from r+ and is bounded from above. Similarly, we also have an interval 
(rrnax " S,rmax + S) with thc property that 



and 



V{r) -uj^ > cA, Vr e (vmax - S, r„iax + S) 
ar 

We repeat the construction of ^11 with the current Q = Qq + AQi and the same 
choice of the functions / and h, which gives 



x: 



(2/' + Ah)\uf + {Ah{V - io') - fV) \uf - ^(/"' + Ah")\uf 

Jr- oo 
[2/Re(u'i/) + (/' + AK) R%{uB)\ , 

with a non-degenerate positive definite integrand on the left hand side. In this case 
however, due to the non super-radiance condition, we may simply subtract a large 
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multiple of the current Qt to handle the boundary terms 



oo r 1 

/ (2/' + Ah)\u'\'' + {AhiV - oj') - fV) \uf - + Ah")\u\' 



+ Buj{uj - uj+)\u\l^^^ + Buj'^\uf\r=, 



Jr+ 



/ [2/ Rc{u'H) + (/' + Ah) Rc{uH) - Blm{uH)] . 



We summarise: 




oo 

< / [-2/ Re{u'H) - (/' + Ah) Rc{uH) + BujIm{uH)] 



11.7. The range. This is a regime in which all frequencies lu, m and Xmi are 
bounded. Recall that in the small a case, a general construction for such frequencies 
was presented in Section 10.2 of this paper. As shown in [18], this construction also 
holds for general a < M in the case m = 0. Moreover, it can be easily adapted 
to cover more generally non-superradiant modes in this regime. To handle the 
remaining bounded superradiant frequencies in this regime, the argument requires 
a modification of Whiting's results on the absence of exponentially growing modes 
in [47]. The details wiU be given in [22]. 

We note that since the azimuthal frequency m can be defined without the Fourier 
transform in time, we can define a priori the projection Pm to the m'th azimuthal 
mode. Clearly, solutions of the wave equation such that Pm'/' = for 1 < |to| < toq 
and large toq are not supported at all in the frequency range C/j,. Thus, for such 
solutions, this appeal to [47] is unnecessary and we have given here a completely 
self-contained discussion of all constructions necessary for the proof. 

11.8. Summation and the local integrated energy decay. We have thus ob- 
tained the inequalities (61) for all frequencies, as required, and we may thus apply 
Section 9.3.3 to obtain 



The precise form of the estimates of Section 11.5 now reveals the nature of the 
physical-space degeneration function x'- The function x must vanish in a neighbor- 
hood of the projection of the 'trapped set' 



Of course, we have really obtained a stronger estimate than (76), namely, we 
obtain the boundedness of the sum of the left hand sides of the versions of (61) 
obtained here, where each microlocalised x degenerates at most at a single point. 
The degeneration of this stronger estimate however, as opposed to (76), can thus 
only be understood microlocally. 

The reader should return to Section 9.3.3 to obtain the remaining statements of 
Theorem 1. 



(76) 
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